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SUMMARY 


Many  empirical  studies  have  found  evidence  that  abrupt  thinning  of  the  crust  can  weaken  the  Lg 
phase,  and  that  contained  sedimentary  basins  or  mountain  ranges  on  the  path  can  completely  or  par¬ 
tially  block  the  phase.  Theoretical  (modeling)  studies  can  provide  insights  for  better  understanding  and 
interpretations  of  the  physical  mechanisms  involved  (Blandford  et  ai,  1992).  This  study  attempts  to 
quantitatively  predict  the  waveguide  effects  using  extensive  two-dimensional  linear  elastic  {P-SV)  finite- 
difference  [LFD]  calculations.  Miscellaneous  waveguide  effects  are  analyzed  and  compared,  both  quali¬ 
tatively  and  quantitatively,  through  visualization  and  spectral  analysis.  A  pure  Lg  wave  packet  is 
injected  into  a  stratified  portion  of  the  grid  as  the  reference  initial  condition  to  trigger  all  LFD  calculations 
for  a  suite  of  heterogeneous  crustal  structures.  An  obvious  advantage  of  this  approach  is  that  the 
effects  due  to  different  types  of  heterogeneity  on  Lg  phase  can  be  isolated  and  evaluated  easily.  Each 
model  contains  a  laterally  heterogeneous  crust  superimposed  on  the  homogeneous  mantle.  The  hetero¬ 
geneities  examined  include  large-scale  lateral  structural  variations  in  the  crustal  waveguide  (such  as 
Moho  uplift,  crustal  thickening,  rugged  free-surface  topography,  embedded  thick  sedimentary  layers, 
etc.),  anelasticity,  and  small-scale  random  heterogeneity. 

The  LFD  results  demonstrate  that  indeed  both  abrupt  changes  in  the  Moho  topography  and  a  thick 
contained  sedimentary  layer  can  cause  Lg  blockage  primarily  through  strong  Lg  -to-Sn  conversions  (with 
actually  different  physical  mechanisms).  Alluvial  basin  structures  cause  strong  reverberations  as  well  as 
very  strong  Lg  -to-Rg  conversions.  If  the  principle  of  seismic  reciprocity  is  invoked,  Rg  -io-Lg  conver¬ 
sion  could  be  an  important  mechanism  of  Lg  excitation  for  certain  structures.  Unlike  the  case  of  Rg 
waves,  a  moderate  free-surface  topography  alone  does  not  seem  to  affect  Lg  propagation  as  much  as 
do  other  types  of  heterogeneity.  Anelasticity  and  small-scale  random  heterogeneity  can  also  explain  the 
Lg  blockage.  An  RMS  velocity  variation  of  8  percent  in  the  whole  crust  is  equivalent  to  a  Qo(Lg )  of 
about  270,  which  would  be  sufficient  to  reduce  the  peak  amplitude  of  1Hz  Lg  waves  by  30  percent  for 
every  100  km  it  traverses.  However,  LFD  calculations  also  indicate  that  these  mechanisms  exhibit  some 
systematic  differences  in  the  couplings  of  Lg  into  other  phases.  Thus,  identifying  the  actual  mechanism 
responsible  for  the  blockage  along  a  specific  path  is  possible. 


Received  for  publication  21  December  1995. 
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WAVEGUIDE  EFFECTS  OF  LARGE-SCALE  STRUCTURAL 
VARIATION,  ANELASTIC  ATTENUATION,  AND 
RANDOM  HETEROGENEITY  ON  SV  Lg  PROPAGATION: 

A  FINITE-DIFFERENCE  MODELING  STUDY 


1.  INTRODUCTION 

The  most  prominent  regional  phase  generated  by  shallow  events  is  the  Lg  phase  identified  by 
Press  and  Ewing  (1952).  The  Lg  phase  receives  considerable  attention  in  recent  years  because  of  its 
potential  use  in  yield  estimation  {e.g.,  Nuttli,  1986)  and  discriminating  between  earthquakes  and  explo¬ 
sions  at  regional  distances  (e.g.,  Blandford.  1981;  Pomeroy  et  a!.,  1982).  Based  on  empirical  observa¬ 
tions,  RMS  Lg  has  been  recognized  as  a  stable  relative-yield  indicator  (Patton,  1988;  Hansen  et  ai, 

1990) .  However,  it  is  also  recognized  that  Lg,  like  Rg,  is  sensitive  to  changes  in  structure  along  its 
path,  which  can  have  deleterious  effects  on  any  role  as  a  discriminant  (see  Lynnes  and  Baumstark, 

1 991 )  or  magnitude  measure  (see  Jih  and  Lynnes,  1 993),  unless  the  propagation  effects  are  accurately 
accounted  for.  Numerous  studies  have  used  the  sensitivity  of  Lg  to  structural  effects  to  map  regions  of 
anomalous  propagation  and  to  try  to  associate  them  with  crustal  structure.  Bias  in  Lg  -based  magnitude 
measurements  were  reported  by  Gregersen  (1984)  in  Greenland.  Very  low  Lg  /S„  amplitude  ratios  have 
been  obsen/ed  after  crossing  the  Tibetan  plateau  (Ruzaikin  et  at.,  1977),  the  North  Sea  grabens  (Gre¬ 
gersen,  1984;  Kennett  et  al. ,  1985),  the  Caspian  Sea,  or  the  Black  Sea  (Levshin  and  Berteussen,  1979; 
Kadinsky-Cade  et  al. ,  1981).  Chinn  et  al.  (1980)  (and  Lynnes  and  BaumstarK  1991)  observed  that  the 
efficiency  of  the  Lg  propagation  is  better  for  paths  parallel  to  the  structural  trend  than  for  paths  in  the 
perpendicular  direction.  Earlier  studies  (Press  aixi  Ewing,  1952;  Oliver  and  Ewing,  1957;  Savarensky 
and  Valdner,  1960)  have  established  that  Lg  does  not  propagate  through  crust  overlain  by  water  deeper 
than  2  km.  On  the  other  hand,  propagation  across  a  marginal  sea  of  continental  shelf  does  not  com¬ 
pletely  quench  Lg ,  but  can  reduce  its  amplitude.  Baumgardt  (1991)  compared  the  crustal  cross-sections 
for  Lg  propagation,  and  he  found  that  the  Lg  blockage  correlates  with  thick  sediment  very  well.  Basi¬ 
cally,  his  obsen/ation  is  that  paths  that  do  not  ctoss  basins  or  for  which  sediments  do  not  vary  by 
greater  than  3  km  exhibit  little  or  no  Lg  blockage  and  scattering.  Baumgardt  (1991)  also  identified  paths 
for  which  the  surface  elevations  and  crustal  thickness  change  substantially,  and  yet  Lg  propagates 
efficiently.  Thus  his  observations  suggest  that  the  near-surface  sediment-thickness  variations  seem  to 
correlate  more  strongly  with  Lg  blockage  than  do  the  crustal-thickness  variations.  Zhang  and  Lay 
(1994)  used  surface  topography  as  a  manifestation  of  the  varying  crustal  structures.  They  found  a 
strong  correlation  between  Sn/Lg  ratios  for  Eurasian  explosions  and  roughness  or  mean  altitude  of  the 
topography  along  the  path,  based  on  a  meager  data  set. 
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Despite  the  long-time  interest,  there  remain  man/  fundamental  questions  about  the  excitation  and 
propagation  of  Lg  waves  phases  to  be  answered.  Numerical  modeling  of  the  Lg  waves  (and  other 
regional  phases)  would  complement  the  empirical  studies  by  providing  more  accurate  interpretations  and 
better  insight  of  the  underlying  physics.  Theoretical  studies  of  Lg  propagation  across  continental  margins 
have  been  conducted  by  Kennett  (1986a),  Maupin  (1989),  Regan  and  Harkrider  (1989),  Cao  and  Muir- 
head  (1993),  and  Gibson  and  Campillo  (1994),  using  different  techniques.  A  simple  geometrical  ray 
theory  can  be  used  to  predict  the  kinematic  property  in  a  qualitative  manner,  as  Kennett  (1986a)  has 
illustrated,  but  would  fail  to  explain  the  dynamic  properties  for  complex  media.  Kennett  (1986b,  1989ab) 
used  a  modal  summation  to  investigate  Lg  propagation  in  stratified  and  weakly  heterogeneous  media 
Mitchell  and  Hwang  (1987)  computed  multi-mode  synthetics  for  1-D  models  with  various  thickness  of 
low-Q  sediments.  Regan  and  Harkrider  (1989)  used  a  hybrid  of  propagator  matrix  and  the  finite-element 
[FE]  methods  to  model  the  SH-Lg  -wave  propagation.  Cao  and  Muirhead  (1993)  applied  a  2-dimensional 
P-SV  finite-difference  method  to  explore  Lg  blockage  and  argue  that  a  water  column  over  the  crust  is  an 
important  factor  in  blocking  Lg  propagation.  Gibson  and  Campillo  (1994)  applied  both  the  dynamic  ray 
tracing  and  the  boundary-integral  equation  methods  to  model  Lg  blockage  in  the  west  Pyrenees  Range, 
near  the  French-Spanish  border.  They  suggest  that  the  unmodeled  scattering  by  small-scale  features 
within  the  lower  crust  is  the  reason  for  the  observed  blockage. 

In  this  study  the  two-dimensional  linear  elastic  (P-SV)  finite-difference  [LFD]  method  is  utilized  to 
model  the  propagation  and  scattering  of  Lg  waves  in  a  suite  of  crustal  models.  Each  model  contains  a 
laterally  heterogeneous  crust  superimposed  on  the  homogeneous  mantle.  The  heterogeneities  exam¬ 
ined  include  anelaslicity,  random  velocity  variation,  rugged  interfaces,  rugged  free-surface  topography, 
as  well  as  embedded  thick  sedimentary  layers.  The  pure  Lg  wave  packet  is  injected  into  a  stratified 
portion  of  the  LFD  grid  as  the  reference  initial  condition  to  trigger  all  LFD  calculations.  An  obvious 
advantage  of  this  approach  is  that  the  effect  due  to  different  types  of  heterogeneity  on  the  Lg  phase  can 
be  easily  isolated  and  evaluated.  A  number  of  different  output  formats  are  available  from  the  LFD  simu¬ 
lation.  In  fact,  one  of  the  advantages  of  the  LFD  method  (and  the  finite-element  method)  is  the  ability  to 
save  particle  displacements  for  any  number  of  receivers  in  the  grid  and  for  all  times.  A  particularly  infor¬ 
mative  way  to  view  wave  propagation  through  the  models  is  by  using  numerical  Schlieren  digrams,  or 
“snapshots”.  Wavefield  snapshots  are  extensively  used  as  a  visualization  aid  throughout  the  study. 

The  procedure  of  implementing  a  pure  Lg  wave  packet  is  described  in  Section  2  and  validated 
with  a  reference  crustal  model.  LFD  experiments  are  then  conducted  with  a  suite  of  canonical  models 
representing  a  variety  of  crustal  structures.  The  geometrical  effects  of  various  structures  on  Lg  propaga¬ 
tion  are  analyzed,  both  qualitatively  and  quantitatively.  In  Sections  4  and  5,  similar  analyses  are  con¬ 
ducted  for  anelastic  and  random  models,  respectively.  All  three  major  categories  of  structures  are  com¬ 
pared  and  equated,  regarding  their  roles  in  causing  Lg  blockage  and  the  coupling  of  Lg  into  other 
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phases.  A  brief  review  on  the  reported  observations  of  Lg  blockage  is  given  in  Section  6.  Rnally,  a 
quantitative  way  of  describing  Lg  blockage  is  explored,  and  an  improved  amplitude/inagnitude  correcting 
procedure  is  presented,  which  includes  a  correction  to  account  for  the  geometric^  block^e. 

2.  LFD  IMPLEMENTATION  OF  Lg  WAVE  PACKET 

The  LFD  method  has  the  advantage  that  the  solution  contains  all  conversions  and  all  orders  of 
multiple  scattering.  It  permits  examinations  of  fairly  general  models  with  arbitrary  complex  variations  in 
material  properties  and  free-surface  geometry.  Furthermore,  it  does  not  require  many  assumptions  com¬ 
monly  invoked  in  other  theoretical  approaches.  The  basic  limitations  to  the  LFD  method  are  the  compu¬ 
tational  cost  and  memory  requirements.  These  constrain  the  size  of  the  grid  and  the  number  of  time 
steps  that  can  be  calculated  in  a  reasonable  time  frame.  Several  non-standard  features  have  been 
incorporated  into  the  LFD  code  used  in  this  study: 

[1]  Explicit  boundary  conditions  for  the  polygonal  free-surface  topography. 

[2]  The  marching-grid  technique  for  stretching  the  propagation  distance  in  the  lateral  direction. 

[3]  A  pure  Lg  wave  packet  suitable  for  LFD  simulations. 

[4J  A  simple  way  to  incorporate  a  causal,  spatially-varying  attenuation  operator  into  LFD  simulations. 

Detailed  descriptions  of  the  explicit  boundary  conditions  arxJ  the  marching  grid  technique  can  be 
found  in  Jih  et  al.  (1988)  and  Jih  (1993a).  Incorporation  of  the  causal,  spatially-varying  attenuation 
operator  will  be  described  in  Section  4  of  this  report.  The  algorithm  generating  the  incident  Lg  wave 
packet  is  analogous  to  the  one  Boore  (1 970)  developed  for  the  fundamental-mode  Love  wave  packet 
Given  a  stratified  crustal  model,  the  Thomson-Haskell  method  is  used  to  compute  the  Lg  dispersion  and 
eigenfunctions  associated  with  this  structure.  The  Lg  eigenfunctions  are  modulated  by  a  Ricker  (1977) 
wavelet  in  the  wavenumber  domain  and  then  converted  to  the  displacement  field  by  the  inverse  Fourier 
transform.  The  resulting  wavefield  is  well  localized  in  both  spatial  and  wavenumber  domains,  and  can  be 
easily  embedded  In  complicated  2-dimensional  models  (for  which  other  forward-modeling  methods  may 
fail)  as  the  initial  condition  to  trigger  the  LFD  simulation.  Ever  since  Boore  (1970)  outlined  the  original 
algorithm  in  his  SH  LFD  calculations,  the  same  procedure  has  been  adopted  frequently  in  elastic  (P-SV) 
LFD  simulations  of  Rg  propagation/scattering  problems.  Just  a  few  studies  along  this  line  include: 
Munasinghe  and  Famell  (1973),  Martel  et  al.  (1977),  Fuyuki  and  Matsumoto  (1980).  Fuyuki  and  Nakano 
(1984),  Levander  (1985),  Toksoz  et  al.  (1986).  McLaughlin  and  Jih  (1986, 1987).  and  Jih  (1993b,  1995). 
The  study  described  herein  is  an  extension  to  Lg  of  the  Rg  analysis  as  presented  In  McLaughlin  and  Jih 
(1986, 1987)  and  Jih  (1993b,  1995).  In  an  earlier  study,  Jih  and  McLaughlin  (1988)  applied  the  principle 
of  reciprocity  to  model  the  excitation  of  Lg  in  various  source  media.  The  goal  of  this  study  Is  to  investi¬ 
gate  various  path  effects  on  Lg ,  not  the  excitation. 
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On  real  seismograms,  the  Lg  phase  often  lacks  a  dear  onset,  but  it  does  have  a  well-defined 
amplitude  maximum  with  a  group  velocity  around  3.5  km/sec.  The  Lg  waves  are  basically  the  interfer¬ 
ence  of  multiply  reflected  S  waves  boundng  back  and  forth  between  the  free  surface  and  the  Moho. 
We  can  envision  that  a  system  of  planar  S  waves  is  set  off  at  the  Moho  at  equal  delay  with  the  same 
post-critical  indination,  while  another  system  of  planar  S  waves  is  set  off  at  the  free  surface  in  a  sym¬ 
metric  manner,  which  corresponds  to  the  bundle  of  reflected  waves.  If  these  two  systems  of  waves  (or 
rays)  are  properly  confined  in  the  same  stratified  region  of  the  crust,  they  will  establish  repetitive 
reflections  (at  both  the  free  surface  and  the  Moho)  with  constructive  interferences.  By  adjusting  the  time 
delay  between  the  consecutive  S  wavefronts,  we  can  obtain  a  denser  (or  coarser)  interference  pattern. 
The  Lg  waves  can  also  be  described  as  the  superposKion  of  many  higher-mode  surface  waves  which 
interfere  to  give  the  complex  observed  waveforms.  For  typical  continental  paths,  the  important  contribu¬ 
tions  come  from  the  stationary  portions  of  the  group  velocity  cun/es  for  about  the  first  10  modes  up  to  1 
Hz  (Knopoff  et  al. ,  1973).  These  modes  correspond  to  waves  trapped  in  the  crustal  wave  guide. 

The  Lg  wave  packet  generated  by  the  procedure  described  earlier  exhibits  all  these  expected 
features.  Rgure  1  gives  the  vertical-component  snapshots  of  Lg  wave  propagation  in  a  single-layer  cru¬ 
stal  model  taken  at  a  temporal  spadng  of  10  seconds.  The  homogeneous  crust  is  30  km  thick  with  P- 
arrd  S-wave  velocities  of  6.2  and  3.58  km/sec,  respectively  (see  Figure  5).  The  Lg  wave  packet  travels 
at  a  group  velocity  of  approximately  3.33  knVsec  (see  Figures  3  and  4),  which  is  appropriate  for  areas 
like  the  western  U.S.  as  suggested  by  Mitchell  and  Hwang  (1987).  Since  no  scattering  mechanism  is 
present,  the  checkerboard-like  pattern  due  to  the  constructive  interference  of  repeatedy  reflected  S 
waves  trapped  in  the  crust  is  retained  at  all  times  (Figure  1).  This  checkerboard-like  pattern  undoubt¬ 
edly  indicates  that  both  interpretations  of  Lg  waves,  either  as  multiply  reflected  S  waves  or  as  higher¬ 
mode  surface  waves,  are  indeed  adequate. 

Like  the  fundamental  mode  Rayleigh  waves  {Rg ),  planar  body  waves  incident  upon  the  free  sur¬ 
face  at  an  oblique  angle  also  exhibit  a  retrograde  rolling.  Since  Lg  waves  are  repetitive  reflections  (at 
both  the  free  surface  and  the  Moho)  of  planar  S  waves,  it  shares  the  same  characteristics.  At  greater 
depths,  however,  the  particle  motion  of  Lg  is  very  different  from  that  of  Rg .  Rgure  2  shows  particle 
motion  plots  of  Lg  synthetics  recorded  at  five  depths  starting  at  the  free  surface  and  2  km  apart.  At  the 
anti-nodes  (for  example,  at  4  km),  the  particle  motion  is  highly  linear. 
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Figure  1.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  a  single-layer  crustal  model  taken 
at  a  temporal  spacing  of  10  seconds.  The  homogeneous  crust  is  30  km  thick  with  P-  and  S-wave  veloci¬ 
ties  of  6.2  and  3.58  km/s,  respectively  (see  Figure  5).  The  Lg  wave  packet  travels  at  a  group  velocity  of 
approximately  3.33  km/s.  Since  no  scattering  mechanism  is  present,  the  checkerboard-like  pattern  due  to 
the  constructive  interference  of  repeatedly  reflected  S  waves  trapped  in  the  crust  is  retained  at  all  times. 


LFD  Simulation  of  Lg  Propagation:  Model  00 
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Figure  2.  Synthetic  seismogram  (right)  and  corresponding  particle  motion  (left)  at  five  depths  for  the  reference 
1  "layer  model  00  (see  Rgure  1).  Vertical  and  horizontal  components  are  shown  in  solid  and  dashed  lines, 
respectively.  The  sensors  are  2  km  apart.  Note  the  retrograde  rolling  at  the  free  surface  (top)  and  the  highly 
polarized  shear  motion  at  antinodes  (for  example,  at  4-km  depth). 
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Figure  3.  Vertical  (top)  and  horizontal  (bottom)  seismic  sections  of  Lg  displacement  synthetics  for  the  reference  Mayer 
model  00.  The  Lg  group  velocity  is  about  3.33  (=  190/[65-8])  km/sec. 
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Rgure  4.  Vertical  (top)  and  horizontal  (txsttom)  seismic  sections  of  Lg  displacement  synthetics  recorded  at  a  depth  of 
10  km  for  the  reference  1 -layer  model  00.  For  the  6th  overtone,  this  depth  is  very  dose  to  an  antinode  and  hence  the 
motion  is  highly  polarized.  In  later  sections  of  this  study,  the  effidency  of  Lg  transmission  is  measured  using  a  deep 
receiver  to  avoid  the  contamination  from  the  Rg  phase. 
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3.  Lg  PROPAGATION  IN  MODELS 

WITH  IRREGULAR  WAVEGUIDE  GEOMETRY 

Numerical  experiments  are  first  conducted  to  evaluate  the  effects  of  variable  crustal  thickness, 
mountainous  relief,  and  sediment  thickness  on  Lg  propagation.  In  order  to  isolate  these  effects,  all  the 
crustal  models  used  in  the  first  set  of  LFD  calculations  are  simple  variations  of  the  single-layer  crustal 
model  {viz.  Model  00).  Table  1  below  gives  a  brief  description  of  these  canonical  models.  The  upper 
portions  of  twenty  five  representative  models  are  shown  in  Rgures  5  through  9. 

Table  1 .  M^or  Structural  Feature(s)  of  Each  Model 


00 

OA 

OAa 

OB 

OBa 

OC 

OCa 

OCa 

OD 

ODa 

ODb 

OE 

OEa 

OEb 

OEc 

OEd 

OF 

OFa 

OFb 

OFc 

2 

2a 

3 
3a 

4 


single-layer  (homogeneous)  crust  of  30  km  thick 

abrupt  Moho  uplift  (30-»6->30) 

gentle  Moho  uplift  (30— >6— >30) 

a  thinning  crust  (extended  Moho  uplift)  (30->6) 

similar  to  OB,  with  a  gentle  transition  (30->6) 

extended  thick  sediment 

similar  to  OC.  with  a  gentle,  linear  transition 

similar  to  OC,  with  a  smooth  transition 

a  contained  sedimentary  basin 

similar  to  OD,  with  a  smooth  boundary 

OD  +  a  sedimentary  layer 

crustal  thickening  (30->45->30,  polygonal  Moho) 

crustal  thickening  (30->45->30,  smooth  Moho) 

crustal  thickening  (30->40->30,  smooth  Moho) 

crustal  thickening  (30->50->30,  smooth  Moho) 

crustal  thickening  (30^50->30,  polygonal  Moho) 

Moho  uplift  +  a  sedimentary  basin 
similar  to  OF,  with  a  smooth  boundary 
Moho  uplift  +  a  sedimentary  basin  of  slow  velocity 
similar  to  OFb,  with  a  smooth  boundary 
a  moderate  topography 

same  topography  as  2  with  isostatic  compensation 
a  moderate  topography 

same  topography  as  3  with  isostatic  compensation 
thickened  crust  with  a  trapezoidal  Moho  and  free  surface 
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Figure  5.  Five  of  the  canonical  crustal  models  tested  in  this  study.  Only  the  uppermost  portion  of 
each  model  is  shown.  The  actual  grid  size  used  in  LFD  calculation  is  1 000  by  500,  corresponding  to 
250  km  by  125  km.  The  material  parameters  of  each  model  are  shown  on  the  right.  In  each  LFD 
experiment,  the  pure  Lg  wave  packet  is  initialized  in  the  stratified  portion  of  the  model,  near  the  left 
edge  of  the  grid.  Model  00  is  the  laterally-invariant  reference  model. 
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Figure  6.  Same  as  Figure  5  except  for  models  OC,  OCa,  OD,  ODa,  and  ODb.  Both  C  and  D 
series  of  models  have  a  thick  sedimentary  layer  in  the  structure.  The  sediment  is  contained  in 
models  of  the  "D”  series. 
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Figure  7.  Same  as  Figure  5  except  for  models  OE,  OEa,  OEb,  OEc,  and  OEd.  All  “E”  series  models 
have  thickened  crust. 
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Figure  8.  Same  as  Figure  5  except  for  models  OF,  OFa,  OFb,  2a,  and  4.  "F”  series  models  are  very 
similar  to  those  of  "D”  series  except  that  a  Moho  uplift  is  added.  Model  2a  has  a  moderate  topogra¬ 
phy  with  a  full  isostatic  compensation.  Model  2  (not  shown)  has  an  identical  tree-surface  topography 
with  a  flat  Moho. 
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Figure  9  shows  the  vertical-component  snapshots  of  Lg  wave  propagation  in  the  model  OA,  which 
has  an  abrupt  Moho  uplift.  The  bulk  of  the  Lg  energy  is  converted  to  Sn  at  the  ascending  crust-mantle 
interface,  as  shown  in  Figure  9  at  20  seconds.  As  a  result,  the  Lg  eigenfunction  is  broken  up  and  only 
a  fraction  of  the  Lg  energy  would  enter  the  oceanic  crust  The  amplitude  of  the  direct  Lg  phase  is 
further  reduced,  due  to  the  geometric  spreading,  when  it  passes  the  shallow  crust  to  enter  the  suddenly 
thickened  continental  crust  (see,  at  40  seconds).  Record  sections  of  surface  synthetics  for  the  model 
OA  are  shown  in  Rgure  10.  While  Rg  remains  intact,  the  Lg  is  effectively  blocked,  although  the  oceanic 
crust  is  only  50  km  long,  the  m^ority  of  Lg  energy  is  converted  to  S„  of  which  some  aosses  the  des¬ 
cending  Moho  and  re-enters  the  crust.  Due  to  the  post-critical  direction,  this  phase  will  stay  in  the  crust 
and  become  Lg  at  a  larger  distance.  This  phase  is  most  prominent  at  43  seconds  on  trace  20,  and  is 
outside  the  expected  “on-time  Lg  ”  window. 

Model  OAa  is  similar  to  OA  except  that  a  gentle,  piecewise-linear  Moho  is  used  instead  (Figure  11). 
The  Lg  -to- Sn  conversion  is  quite  clear  at  20  and  30  seconds.  The  energy  which  crosses  the  oceanic 
crust  as  the  3rd  overtone  (see,  at  30  seconds)  does  not  touch  the  free  surface  as  the  crust  thickens. 
Instead,  it  propagates  along  the  descending  Moho  (see,  at  40  and  50  seconds)  before  it  reaches  the  fiat 
Moho  and  starts  to  bounce  back  to  the  ground  (see,  at  60  seconds).  The  “Lg  blockage”  occurs  because 
only  S„  and  Rg  could  reach  ground  sensors  in  the  oceanic-to-continental  transition  zone.  Rgure  12 
shows  that  the  dominating  phase  of  this  model  is  the  “early  Lg  ”  with  a  group  velocity  about  4.0  km/sec, 
followed  by  Rg  waves  with  a  group  velocity  of  3.0  km/sec.  There  is  no  prominent  phase  arriving  at  the 
expected  “on-time”  Lg  window  of  3.3  km/sec. 

Model  OB,  shown  in  Figure  13,  has  a  crust  that  thins  abruptly  from  a  thickness  of  30  km  to  6  km. 
As  it  does  in  the  “A”-series  of  models,  a  significant  fraction  of  Lg  energy  enters  the  deeper  mantle  and 
detaches  from  the  Moho  (see,  30,  40,  and  50  seconds  of  Rgure  13).  This  energy  never  reaches  the 
surface  sensors,  and  as  a  result,  the  most  prominent  phase  is  the  stabilized  3rd  overtone,  which  stays  in 
the  thin  crust.  Record  sections  of  this  model  (Rgure  14)  indicate  that  no  “Lg  blockage”  would  be 
observed  here  since  the  3rd  overtone  has  a  peak  amplitude  and  a  group  velocity  very  similar  to  those  of 
the  “on-time  Lg  ”  in  the  reference  model.  Changing  the  Moho  shape  does  not  affect  these  observations 
(see  Rgures  15  and  16). 

All  “C”'Series  models  have  a  sedimentary  layer  16  km  thick.  The  sediment  slows  down  the  Lg 
waves,  as  expected  (Rgures  17  and  19).  However,  this  structure  does  not  generate  a  Sn  conversion 
as  strong  as  in  other  models  with  an  irregular  Moho  or  contained  basin  (Rgures  18  and  20). 

For  models  of  the  “D”  series,  because  of  the  contained  structure  of  the  basin,  the  Lg  wave  that 
tunnels  underneath  the  basin  eventually  recovers  characteristics  like  a  multiply  bouncing  SmS  (Rgure 
21).  This  phase  has  the  same  group  velocity  of  the  “on-time  Lg  ”  on  those  surface  sensors  co-located 
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at  the  bounce  points.  The  amplitude  is  smaller,  however.  A  more  prominent  phase  comes  from  the 
energy  that  first  enters  the  basin  and  then  leaks  out  from  the  ascending  edge  of  the  basin.  The  comer 
where  the  pinched  basin  ends  also  radiates  some  body  waves.  Rgure  22  illustrates  that  the  contained 
basin  structure  removes  some  energy  of  the  “on-time  Lg  ”  and  converts  it  to  the  so-called  “late  Lg  ” 
which  has  the  same  phase  velocity  as  the  reference  Lg  but  is  significantly  delayed.  Depending  on  the 
location  of  surface  sensors,  some  seismographs  will  observe  either  complete  or  partial  blockage  of  the 
“on-time  Lg  ”.  Rgures  23-26  show  results  for  similar  structures  (models  ODa  and  ODb). 

Model  OE  has  a  thickened  crust  and  a  very  strong  Lg  -to-S„  conversion  occurs  along  the  ascend¬ 
ing  interface  (see,  30,  40,  and  50  seconds  of  Rgure  27).  Horizontal-component  snapshots  (not  shown) 
indicate  that  some  focusing  of  energy  occurs  at  the  comer  where  the  Moho  turns  flat.  This  comer 
becomes  a  secondary  point  source  radiating  upgoing  body  waves.  Sensors  right  above  this  comer  would 
detect  a  prominent  arrival  due  to  a  nearly  vertical  incidence  of  shear  waves  (see  traces  16  and  17  of 
Figure  28).  Rgures  29-36  show  variations  of  thickened  sections  of  crust  (models  OEa  through  OEd). 

“F’-series  models  have  both  a  sedimentary  basin  and  a  Moho  uplift.  Rgure  37  shows  the  basic 
geometry.  The  Moho  uplift  causes  a  strong  S„  conversion  (see,  at  20  seconds  of  Rgure  38).  The 
pinched  end  of  the  basin  radiates  body  waves  outward,  and  some  of  those  rays  leak  into  deeper  mantle. 
Figure  39  shows  the  horizontal-component  snapshots  of  Lg  wave  propagation  in  the  model  OFa  which  is 
similar  to  OF  except  for  a  smooth  interfaces.  Observations  similar  to  those  with  the  model  OF  can  be 
made.  The  pinched  end  of  the  basin  radiates  body  waves  outward,  and  some  of  those  rays  leak  into 
deeper  mantle  (see,  at  50  seconds).  Some  direct  P  waves  are  also  present  (see,  at  40  seconds).  Direct 
P  waves  radiated  from  the  basin’s  terminating  edge  can  be  seen  on  the  horizontal  components  of  Rgure 
40.  Figure  41  -44  show  the  geometry  and  results  for  some  other  variations  (models  OFb  and  OFc). 

Model  2  has  a  flat  Moho  and  a  moderate  free-surface  topography  (Figure  45).  The  Lg  -to-S„ 
conversion  is  relatively  weak  (see  Figure  46).  This  is  very  different  from  the  case  of  Rg  where  scatter¬ 
ing  by  rough  topography  would  cause  very  strong  apparent  attenuation.  Model  2a  has  a  fully  compen¬ 
sated  Moho  added  to  the  model  (Rgure  47).  The  fully  compensated  Moho  gives  results  similar  to  those 
from  the  flat  Moho  of  model  2.  Figures  47  through  54  give  results  from  models  2a,  3,  3a,  3a,  and  4, 
which  are  all  similar  to  model  2. 

Figures  55  and  56  display  the  vertical-component  synthetic  seismograms  recorded  at  3  depths 
(0km,  10km,  and  40km)  and  started  eit  20  seconds  for  models  00,  OA,  OAa,  OB,  OBa,  OC,  OCa,  OD,  ODa, 
and  ODb.  Rgures  57  and  58  show  the  same  for  models  OE,  OF,  OFa,  OFc,  and  2a  starting  at  time  zero. 
Note  that  although  the  “on-time  Lg  ’’  on  surface  synthetics  of  models  OB  and  OBa  (bottom  left)  have  a 
peak  amplitude  in  the  same  order  as  that  of  the  reference  model,  they  are  not  the  same  overtone  of  the 
Rayleigh  mode.  The  moderate  surface  and  Moho  topography  in  the  model  2a  only  affect  Lg  waveform 
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characteristics  mildly.  Note  that  model  OFc  excites  very  strong  Rg  phase,  very  similar  to  the  model  ODb. 
Among  all  structures  tested,  “E”  series  of  models  (with  a  thickened  crust)  excites  the  strongest  S„ 
conversion  (see  Figure  27). 

It  is  interesting  to  note  that  all  models  of  the  “F’  series  exhibit  significant  Lg  -lo-Pg  conversion. 
The  peak  horizontal  amplitude  on  trace  No.25  ranges  from  10%  to  15%  of  that  of  the  original  Lg  wave 
train.  To  a  lesser  extent,  models  of  “C”  and  “D”  series  as  well  as  models  OA  and  OBa  also  generate 
the  same  conversion. 

For  each  model,  the  synthetic  seismograms  are  recorded  at  three  linear  arrays  at  three  different 
depths:  0  km,  10  km,  and  40  km.  Several  parameters  are  measured  to  infer  the  relationship  between  the 
major  geological/geophysical/geometrical  features  and  possible  blockage  and  coupling.  The  ampli¬ 
tude  spectrum  of  a  surface  synthetic  seismogram  far  away  from  the  heterogeneous  structure  is  com¬ 
puted  using  a  window  of  12.8  seconds  (512  data  points)  around  the  group  velocity  of  3.33  km/sec.  Using 
the  spectrum  of  model  00  as  a  reference,  the  amplitude  ratio  at  1  Hz  is  denoted  as  “BK”  in  Table  2  to 
quantify  the  Lg  blockage.  This  ratio  is  also  used  to  compute  the  equivalent  apparent  attenuation  of  the 
“on-time  Lg".  Both  the  quality  factor  Qq  and  the  attenuation  coefficient  y  are  printed  in  Table  2.  To 
measure  the  total  transmitted  Lg  energy,  a  sensor  at  10  km  depth  is  used.  The  spectral  amplitude  ratio 
at  1  Hz  computed  with  this  sensor  is  denoted  as  ‘T”  in  Table  2.  The  converted  S„  energy  is  measured 
with  a  sensor  at  a  depth  of  40  km,  below  the  flat  Moho.  Another  sensor,  which  is  near  the  bottom  of  the 
model  and  under  the  major  heterogeneous  feature,  records  the  waveform  representing  the  converted 
teieseismic  energy.  Except  for  the  quantity  BK,  all  other  ratios  are  normalized  to  the  lOkm-deep  syn¬ 
thetic  of  the  model  00  as  the  reference.  The  results  in  Table  2  are  sorted  by  the  the  major  structural 
features  of  the  models.  There  appear  to  be  some  systematic  patterns  that  are  very  interesting. 
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Table  2.  Lg  Propagation  Statistics  of  Deterministic  Models 


Model 

BK 

T 

SniLg 

far-P/Lg 

far-S/Lg 

Qo 

y(1  Hz) 

OA 

0.349 

0.346 

0.164 

0.110 

0.105 

91 

0.0105 

OAa 

0.416 

0.689 

0.216 

0.083 

0.087 

109 

0.0088 

OAb 

0.384 

0.624 

0.233 

0.080 

0.079 

100 

0.0096 

OB 

1.035 

0.546 

0.109 

0.103 

OBa 

0.848 

0.482 

0.041 

0.044 

581 

0.0016 

OC 

0.354 

0.790 

0.084 

0.061 

0.048 

92 

0.0104 

OCa 

0.237 

0.521 

0.076 

0.026 

0.010 

66 

0.0144 

OCb 

0.307 

0.496 

0.065 

0.051 

0.034 

81 

0.0118 

OD 

0.595 

0.565 

0.161 

0.075 

0.084 

185 

0.0052 

ODa 

0.317 

0.577 

0.185 

0.073 

0.096 

84 

0.0115 

ODb 

0.508 

0.409 

0.102 

0.083 

0.086 

142 

0.0068 

OE 

0.977 

0.779 

0.468 

0.046 

0.051 

4136 

0.0002 

OEa 

0.963 

0.723 

0.459 

0.038 

0.082 

2534 

0.0004 

OEb 

0.967 

0.649 

0.201 

0.032 

0.050 

2821 

0.0003 

OEc 

0.983 

0.719 

0.478 

0.052 

0.070 

5558 

0.0002 

OEd 

0.923 

0.741 

0.318 

0.058 

0.039 

1191 

0.0008 

OF 

0.228 

0.580 

0.177 

0.058 

0.062 

65 

0.0148 

OFa 

0.789 

0.595 

0.225 

0.076 

0.084 

405 

0.0024 

OFb 

0.318 

0.533 

0.157 

0.070 

0.087 

84 

0.0115 

OFc 

0.473 

0.560 

0.194 

0.078 

0.098 

128 

0.0075 

2 

1.023 

0.574 

0.087 

0.058 

0.127 

_ 

2a 

1.020 

0.717 

0.152 

0.059 

0.138 

_ 

_ 

3 

0.875 

0.829 

0.091 

0.061 

0.145 

721 

0.0013 

3a 

1.111 

0.715 

0.133 

0.070 

0.150 

— 

— 

4 

1.738 

0.408 

0.319 

0.050 

0.072 

— 

— 

BK;  sp6Ctral  amplitudG  ratio  at  1  Hz  of  “on-tim©  Lg  ”  (3S3  kin/sGc)  to  that  of  the  reference  model. 

T;  spectral  amplitude  ratio  at  1  Hz  of  transmitted  Lg  (recaded  at  a  depth  of  10  km)  to  that  of  the  reference  model. 

Sn/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  S„  to  the  reference  Lg  . 

far-P/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  teleseismic  P  to  the  reference  Lg  . 

far-S/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  teleseismc  S  to  the  reference  Lg  . 


For  seismic  stations  in  northern  Scandinavia,  the  1  Hz  Lg  waves  from  Novaya  Zemlya  are 
severely  blocked.  Based  on  the  coda  Q  method,  Nuttli  (1988)  deduced  a  Qq  value  of  252  for  the  path 
from  Novaya  Zemlya  to  WWSSN  station  KEV  (Kevo,  69.755»N,  27.007“E).  Using  a  joint  inversion 
method  with  Lg  amplitude  measurements.  Jih  et  al.  (1995)  point  out  that  this  Qq  value  appears  to  be 
appropriate.  Assuming  that  the  coda  Q  method  really  provides  an  independent  and  consistent  measure 
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of  Lg  Q,  then  there  should  be  a  Lg  -mb  bias  of  0.26  magnitude  unit  [m.u.]  between  the  Eastern  U.S. 
and  Novaya  Zemlya.  The  Qq  value  of  252  for  station  KEV  actually  includes  the  gross  reduction  effects 
on  Lg  amplitude  due  to  both  the  intrinsic  attenuation  and  the  geometrical  blockage,  of  which  the  latter 
was  not  fully  studied  or  documented  at  the  time  Nuttli  estimated  the  yield  for  Novaya  Zemlya  explosions. 
The  continental  path  crossing  Barents  Shelf  is  a  classical  example  of  “Lg  blockage”  (Baumgardt,  1991). 
For  blockage  of  the  “on-time  Lg  ”  to  be  noticeable  at  similar  distance,  it  is  reasonable  to  expect  to  have 
a  Qo  value  no  more  than  270,  or  equivalently,  a  BK  no  greater  than  approximately  0.7. 

Using  a  BK  of  0.7  as  the  threshold,  we  see  that  all  models  of  the  “A”,  “C”,  and  “D”  series  as  well 
as  most  of  the  “F”  series  would  block  Lg ,  while  none  of  the  “B”  and  “E”  models  or  any  model  with  a 
mild  topography  would  block  Lg  .  All  “E”-type  of  models,  which  have  a  thickened  crust,  exhibit  very  little 
attenuation.  Most  of  the  topographic  models  (that  is,  models  2,  2a,  3a,  and  4)  actually  exhibit  an 
amplification,  instead  of  reduction,  in  Lg  amplitude.  The  thinning  of  the  crust  is  an  important  mechanism 
to  block  Lg  (Rgures  8  and  10),  as  reported  in  many  observational  studies.  However,  the  thinning  itself 
has  to  be  of  finite  span  in  the  lateral  direction.  The  two  “B”  series  of  models  either  show  an 
amplification  or  have  a  BK  larger  than  0.85  (see  Figures  12  and  14). 

The  comparison  of  Lg  coupling  to  teleseismic  phases  is  also  described  in  Table  2.  In  terms  of  Lg  - 
\o-S„  conversion,  models  OB  and  OBa  are  the  most  efficient  ones,  followed  by  four  models  with  thick¬ 
ened  crust,  OE,  OEa.  OEc,  and  OEd.  Another  “E”-type  model,  OEb,  with  an  elliptical  Moho  shape  as  in 
model  OEa  but  only  40  km  thick,  gives  a  much  weaker  S„  excitation.  On  the  other  hand,  all  these  “E” 
series  models  have  very  comparable  efficiency  in  Lg  transmission  -  either  measured  at  the  free  surface 
or  at  a  depth  of  10  km.  Model  4  has  a  trapezoidal  mountain  (5  km  high)  on  the  top  and  a  trapezoidal 
intrusion  of  10  km  in  the  bottom  of  the  crust  Thus  its  overall  increase  in  thickness  is  identical  to  that  of 
models  OE  and  OEa.  While  the  level  of  S„  conversion  is  in  the  same  range  as  other  “E”  models,  the  Lg 
transmission  is  actually  amplified.  If  the  models  are  sorted  by  the  corresponding  levels  of  S„  conver¬ 
sion,  we  have  {0B,0Ba}  >  {0E,0Ea.0Ec,0Ed}  »  4  >  OA,OEb  >  OF*  >  0D*,2a,3a  >  2,3  >  OC*. 
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LFD  Simulation  of  Lg  Propagation:  Model  OA 

Figure  9.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  the  model  OA,  which  has  an 
abrupt  Moho  uplift.  The  bulk  of  Lg  energy  is  converted  to  S„  at  the  ascending  crust-mantle  interface 
(see.  at  20  seconds).  As  a  result,  the  Lg  eigenfunction  is  broken  up  and  only  a  fraction  of  the  Lg  energy 
would  enter  the  oceanic  aust.  The  amplitude  of  the  direct  Lg  phase  is  further  reduced  (due  to  the 
geometric  spreading)  when  it  passes  the  shallow  crust  to  enter  the  suddenly  thickened  continental  crust  (at 
40  seconds).  Although  the  oceanic  crust  is  only  50  km  long,  it  effectively  blocks  the  Lg . 
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Figure  10.  Record  sections  of  surface  synthetics  for  the  model  OA.  While  Rg  is  intact,  the  Lg  is  effectively  blocked. 
The  majority  of  Lg  energy  is  converted  to  S„  of  which  some  crosses  the  descending  Moho  and  re-enters  the  crust.  Due 
to  the  post-critical  direction,  this  phase  will  stay  in  the  crust  and  become  Lg  sA  a  larger  distance.  This  phase  is  most 
prominent  at  43  seconds  on  trace  20.  and  is  outside  the  expected  “on-time  Lg  ”  window. 


22 


Rlstep  0 

VEFTT  0.000SOC,  1  800  1  500 

XTicMarfcs  10  km;  Z  Tic  Mark  =  10km 
-O^E-fOO,  0^6E+00 

FDatop  400 

VERT10.000sec,  1  800  1  500 
XTIcMai1c=  10  km;  2  Tic  Marks  10km 
-0^7E+00,  O^E-^00 

FDstep  800 

VERT  20X)00s«C,  1  800  1  500 

XTlcMarks  10km;  ZTIcMarks  10km 
-0^E400,  O^E+OO 

FDatep  1200 

VERT80.000sec,  1  800  1  500 
XTIcMailcs  10  km;  Z  Tic  Marks  lOkm 
-02BE*00,  O^E+00 

FDstap  1600 

VERT  AOJOOOsbc,  1  800  1  500 
XTIcMarks  10km;  ZTIcMarks  10km 
-OJ26E*O0,  0J22EM 

FDstop  2000 

VERT50.000MC,  1  800  1  500 

XTIcMarks  10km;  ZTIcMarks  10km 
-0.19E+00,  0.18E44M) 

FDstep  2400 

VERT  60iX)0aec,  1  800  1  500 
XTIcMarks  10km;  ZTIcMarks  10km 
^^E+00,  0J22E*O0 


LFD  Simulation  of  Lg  Propagation:  Model  OAa 

Figure  11.  Same  as  Figure  9  except  for  the  model  OAa.  which  has  a  gentle  continental-to-oceanic  transi¬ 
tion.  The  Lg  -to-Sn  conversion  is  quite  clear  at  20  and  30  seconds.  The  energy  that  crosses  the  oceanic 
crust  as  the  3rd  overtone  (see,  at  30  seconds)  does  not  touch  the  free  surface  as  the  crust  thickens. 
Instead,  it  propagates  along  the  descending  Moho  (at  40  and  50  seconds)  before  it  reaches  the  flat  Moho 
and  starts  to  bounce  back  to  the  ground  (see,  at  60  seconds).  The  "Lg  blockage”  occurs  because  only 
S„  and  Rg  could  reach  ground  sensors  in  the  oceanic-to-continental  transition  zone. 
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Synthetics  at  0km  Depth,  Model  OAa 
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Rgure  12.  Seismic  sections  of  surface  synthetics  for  the  model  OAa.  The  dominating  phase  is  the  “early  4  ”  with  a 
group  velocity  about  4.0  km/sec,  followed  by  Rg  waves  with  a  group  velocity  of  3.0  km/sec.  There  is  no  prominent 
phase  arriving  at  the  expected  “on-time”  4  window  at  3.3  km/sec. 
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LFD  Simulation  of  Lg  Propagation:  Model  OB 

Figure  13.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OB  in  which  the  crustal 
thickness  decreases  abruptly  from  30  to  6  km.  A  significant  fraction  of  Lg  energy  enters  the  deeper  man¬ 
tle  and  detaches  from  the  Moho  (see  30,  40.  and  50  seconds).  This  energy  never  reaches  the  surface 
sensors,  and  as  a  result,  the  most  prominent  phase  is  the  stabilized  3rd  overtone,  which  stays  in  the  thin 
crust. 
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Rgure  14.  Seismic  sections  of  surface  synthetics  for  model  OB.  The  most  prominent  phase  is  the  3rd  overtone  con¬ 
verted  from  the  original  Lg .  Note  that  no  “Lg  blockage”  would  be  observed  here  since  this  3rd  overtone  has  a  peak 
amplitude  and  a  group  velocity  very  similar  to  those  of  the  "on-time  Lg  "  in  the  reference  model  (see  Figure  3). 
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LFD  Simulation  of  Lg  Propagation:  IVIodel  OBa 

Figure  15.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OBa  with  a  crustal  thick¬ 
ness  gradually  decreasing  from  30  to  6  km.  As  in  model  OB,  the  Lg  is  only  partially  blocked  with  some 
energy  traveling  within  the  thin  crust  as  the  3rd  overtone.  Much  of  the  energy  leaks  into  deeper  mantie 
apparently  in  the  form  of  body  waves  radiating  from  the  corner  where  the  thin  crust  starts,  which  is  easier 
to  see  on  the  horizontal  component. 
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Figure  16.  Seismic  sections  of  surface  synthetics  for  model  OBa.  The  most  prominent  phase  on  trace  No.  26  is  the  3rd 
overtone  converted  from  the  original  Lg .  which  arrives  “on-time”  and  hence  only  a  partial  blockage  is  observed.  A 
strong  Lg  -to-P  conversion  at  the  comer  generates  a  very  prominent  arrival  on  the  horizontal  component.  This  phase  dies 
out  very  fast  on  the  surface  recordings. 
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LFD  Simulation  of  Lg  Propagation:  Nlodei  OC 

Figure  17.  The  verticai-component  snapshots  of  Lg  wave  propagation  in  model  OC.  which  has  a  sedi¬ 
mentary  layer  16  km  thick.  The  sediment  siows  down  the  Lg  waves,  as  expected.  However,  this  structure 
does  not  generate  a  S„  conversion  as  strong  as  in  other  models  with  irregular  Moho  or  contained  basin. 
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Rgure  18.  Seismic  sections  of  surface  synthetics  for  model  OC. 
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Figure  19.  Same  as  Figure  17  except  that  the  basin  has  a  gentle  edge. 
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Figure  20.  Seismic  sections  of  surface  synthetics  for  model  OCa. 
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LFD  Simulation  of  Lg  Propagation:  Model  OD 

Figure  21.  Because  of  the  contained  structure  of  the  basin,  the  Lg  that  tunnels  underneath  the  basin 
eventually  recovers  characteristics  like  a  multiply  bouncing  SmS.  This  phase  has  the  same  group  velocity 
as  the  “on-time  Lg  ”  on  those  surface  sensors  co-located  at  the  bounce  points.  The  amplitude  is  smaller, 
however.  A  more  prominent  phase  comes  from  the  energy  that  first  enters  the  basin  and  then  leaks  out 
from  the  ascending  edge  of  the  basin.  The  comer  where  the  pinched  basin  ends  also  radiates  some  body 
waves. 
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Rgure  22.  The  contained  basin  structure  removes  some  energy  from  the  “on-time  Lg  ”  and  converts  it  to  the  so-called 
“late  Lg  ”  which  has  the  same  phase  velocity  as  the  reference  Lg  but  is  significantly  delayed.  Depending  on  the  location 
of  surface  sensors,  some  sensors  will  observe  either  complete  or  partial  blockage  of  the  "on-time  Lg 
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LFD  Simulation  of  Lg  Propagation:  Model  ODa 

Figure  23.  Same  as  Figure  21  except  for  an  elliptical  basin.  The  results  are  very  similar  to  those  of 
model  OD.  Thus  the  details  of  the  basin  shape  appear  to  be  not  that  important. 
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Synthetics  at  0km  Depth,  Model  ODa 

HORl,  same  scale,  decimated  by  5 

Figure  24.  Similar  to  Figure  22  except  for  an  elliptical  basin.  Although  the  results  are  very  similar  to  those  of  model  OD, 
the  blockage  of  the  “on-time  Lg  ”  is  more  apparent.  Also  note  the  strong  backscattering  at  the  terminating  edge  of  the 
basin. 
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LFD  Simulation  of  Lg  Propagation:  Model  ODb 

Figure  25.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  the  basin  model  ODb.  As  in 
models  OD  and  ODa,  the  bodywave  coupling  occurs  at  the  pinched  end  of  the  basin  where  part  of  the 
energy  that  traverses  the  basin  starts  to  leak  out  (at  40  and  50  seconds).  This  suggests  that  even  in  a 
model  with  no  lateral  variation  in  the  Moho,  the  lateral  structural  heterogeneity  in  the  uppermost  crust  still 
can  produce  mantle  phases  through  scattering.  The  pinched  edge  is  essential  in  this  mechanism. 
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Figure  26.  Similar  to  Rgure  22  except  that  a  thin,  slow  sedimentary  layer  is  added.  Both  the  Rg  phase  and  the  "late 
Lg  ”  become  prominent  in  this  case.  The  "on-time  Lg  "  wave  is  weakened  and  only  certain  surface  sensors  at  the  right 
locations  can  observe  it. 
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LFD  Simulation  of  Lg  Propagation:  Model  OE 

Figure  27.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OE  with  a  thickened  crust. 
A  very  strong  Lg  -to-S„  conversion  occurs  along  the  ascending  interface  (see.  30,  40,  and  50  seconds). 
Horizontal-component  snapshots  (not  shown)  indicate  that  some  focusing  of  energy  occurs  at  the  corner 
where  the  Moho  turns  flat.  This  corner  becomes  a  secondary  point  source,  radiating  upgoing  body  waves. 
Sensors  right  above  this  corner  would  detect  a  prominent  arrival  due  to  a  nearly  vertical  incidence  of  shear 
waves  (see  traces  No.  1 6  and  1 7  of  Figure  28). 
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Rgure  28.  Seismic  sections  of  model  OE,  which  has  a  thickened  crust  The  corner  where  the  Moho  turns  flat  is  a 
secondary  point  source  radiating  upgoing  body  waves.  The  prominent  arrival  on  traces  16  and  17  has  a  very  large 
apparent  velocity,  due  to  a  nearly  normal  incidence.  The  “on-time  Lg  "  is  only  partially  blocked. 
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Figure  29.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OEa  with  a  thickened 
crust. 


LFD  Simulation  of  Lg  Propagation:  Model  OEa 
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Rgure  30.  Seismic  sections  of  model  OEa,  which  has  a  thickened  cmst.  The  comer  where  the  Moho  turns  fiat  is  a 
secondary  point  source  radiating  upgoing  body  waves.  The  prominent  arrival  on  traces  16  and  17  has  a  very  large 
apparent  velocity,  due  to  a  nearly  normal  incidence.  The  “on-time  Lg  ”  is  only  partially  blocked. 
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LFD  Simulation  of  Lg  Propagation:  Model  OEb 

Figure  31.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OEb.  with  a  thickened 
crust. 
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Figure  32.  Seismic  sections  cf  model  OEb,  which  has  a  thickened  crust  The  comer  where  the  Moho  turns  flat  is  a 
secondary  point  source  radiating  upgoing  body  waves.  The  prominent  arrival  on  traces  16  and  17  has  a  very  large 
apparent  velocity,  due  to  a  nearly  normal  incidence.  The  “on-time  Lg  ”  is  only  partially  blocked. 
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LFD  Simulation  of  Lg  Propagation:  Model  OEc 

Figure  33.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OEc.  with  a  thickened 
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LFD  Simulation  of  Lg  Propagation:  Model  OEd 

Figure  35.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OEd,  with  a  thickened 
crust 
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Rgure  36.  Seismic  sections  of  model  OEd,  which  has  a  thickened  crust  The  comer  where  the  Moho  turns  flat  is  a 
secondary  point  source  radiating  upgoing  body  waves.  The  prominent  arrival  on  traces  16  and  17  has  a  very  large 
apparent  velocity,  due  to  a  nearly  normal  incidence.  The  “on-time  Lg  "  is  only  partially  blocked. 
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LFD  Simulation  of  Lg  Propagation:  Model  OF 

Figure  37,  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OF.  which  has  both  a 
sedimentary  basin  and  a  Moho  uplift.  The  Moho  uplift  causes  a  strong  Sn  conversion  (see,  at  20 
seconds).  The  pinched  end  of  the  basin  radiates  body  waves  outward,  and  some  of  those  rays  leak  into 
deeper  mantle. 
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Figure  38.  Seismic  sections  of  model  OF,  which  has  a  Moho  uplift  and  a  sedimentary  basin. 
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LFD  Simulation  of  Lg  Propagation:  Model  OFa 

Figure  39.  The  horizontal-component  snapshots  of  Lg  wave  propagation  in  model  OFa.  which  has  both  a 
sedimentary  basin  and  a  Moho  uplift.  Observations  similar  to  those  with  the  model  OF  can  be  made.  The 
pinched  end  of  the  basin  radiates  body  waves  outward,  and  some  of  those  rays  leak  into  deeper  mantle 
(see.  at  50  seconds).  Some  direct  P  waves  are  also  present  (see.  at  40  seconds). 
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Synthetics  at  Okm  Depth,  Model  OFa 

HORI,  same  scale,  dec^ated  by  5 


Rgure  40.  Seismic  sections  of  model  OFa.  which  has  a  Moho  uplift  and  a  sedimentary  basin.  Partial  blockage  of  "on- 
time  Lg  ”  is  observable.  Much  of  the  energy  is  delayed  because  of  the  sedimentary  basin.  Direct  P  waves  radiated  from 
the  basin’s  terminating  edge  can  be  seen  on  the  horizontal  components. 
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LFD  Simulation  of  Lg  Propagation:  Model  OFb 

Figure  41.  The  horizontal-component  snapshots  of  Lg  wave  propagation  in  model  OFb,  which  has  two 
sedimentary  layers  and  a  Moho  uplift.  Observations  similar  to  those  with  the  models  OF  and  OFa  can  be 
made:  the  basin  edge  acts  like  a  secondary  point  source.  A  prominent  SmS-mode  of  propagation  is  quite 
clear  at  50  seconds.  The  visibility  of  this  phase  is  distance-dependent. 
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Rgure  42.  Seismic  sections  of  model  OFb,  which  has  a  Moho  uplift  and  two  sedimentary  layers.  “On-time  Lg  "  is  essen¬ 
tially  blocked  whereas  two  groups  of  “late  Lg  ”  waves  as  well  as  the  further  delayed  Rg  are  present.  The  prominent 
phase  arriving  between  55  and  70  seconds  on  traces  20  through  24  is  a  SmS  phase  originated  from  the  basin  edge. 
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LFD  Simulation  of  Lg  Propagation:  Model  OFc 

Figure  43.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  model  OFc.  which  has  two  sedi¬ 
mentary  layers  and  a  Moho  uplift  Observations  similar  to  those  with  model  OFb  can  be  made:  the  basin 
edge  acts  like  a  secondary  point  source.  A  prominent  SmS-mode  of  propagation  is  quite  clear  at  50  and 
60  seconds.  The  visibility  of  this  phase  is  distance-dependent,  however. 
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Figure  44.  Seismic  sections  of  model  OFc,  which  has  a  Moho  uplift  and  two  sedimentary  layers.  The  most  prominent 
phase  before  Rg  is  the  SmS  radiated  from  the  basin  edge,  which  can  be  seen  at  certain  receivers,  depending  on  the 
distance  (see  traces  No.  21  through  26). 
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LFD  Simulation  of  Lg  Propagation:  Model  2 


Figure  45.  The  vertical-component  snapshots  of  Lg  wave  propagation  of  model  2,  which  has  a  flat  Moho 
and  a  moderate  free-surface  topography.  The  Lg  -to-S„  conversion  is  relatively  weak.  This  is  very 
different  from  the  case  of  Rg  where  scattering  by  rough  topography  would  cause  very  strong  apparent 
attenuation. 
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Figure  46.  Seismic  sections  of  model  2.  which  has  a  moderate  free-surface  topography.  The  scattering  effects  of 
topography  on  Lg  are  not  as  strong  as  those  due  to  large-scale  structural  discontinuities. 
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LFD  Simulation  of  Lg  Propagation:  Model  2a 


Figure  47.  The  vertical-component  snapshots  of  Lg  wave  propagation  of  model  2a,  in  which  the  Moho  is 
fully  compensated.  The  elevation  change  of  this  topographic  profile  is  relatively  small  compared  to  the 
thickness  of  the  homogeneous  crust,  and  hence  the  transmission  of  Lg  is  still  fairly  efficient. 
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Figure  48.  Seismic  sections  of  model  2a,  which  has  a  moderate  free-surface  topography  and  a  fully  compensated 
Moho.  The  scattering  effects  of  topography  on  Lg  are  not  as  strong  as  those  due  to  large-scale  structural  discontinuities. 


60 


FDstap  0 

VERT  .OOOsec,  1  800  1  500 
XTicMarkr  10  km;  Z  Tic  Mark  =  10km 
^^6E+00,  0^6E4^ 


FDatop  400 

VERT10.0008ac,  1  800  1  500 

XTIcMark=  10km;  ZTlcMarks  10km 
-O^Ei>00,  O^E+00 


FOstep  800 

VERT  20.0008ec,  1  800  1  500 

XTIcMarks  10  km;  Z  Tic  Mark  =  10km 
^^E+00,  O^E-i^ 


FDstap  1200 

VERT30.000sac,  1  800  1  500 
XTIcMarfcs  10  km;  Z  Tic  Mark  =  10km 
-O^Ei^,  0J3OB-M 


FDstap  1600 

VERT40.000sac,  1  800  1  500 
XTIcMarkr  10  km;  Z  Tic  Marks  10km 
-O^E-i^,  O^E^ 


FDstap  2000 

VERTSOUXXteac,  1  800  1  500 
XTIcMarfcs  10  km;  Z  Tic  Mark  =  10km 
-O^IE-^00,  0^1  E^ 


LFD  Simulation  of  Lg  Propagation:  Model  3 


Figure  49.  The  vertical-component  snapshots  of  Lg  wave  propagation  of  model  3.  The  elevation  change 
of  this  topographic  profile  is  relatively  small  compared  to  the  thickness  of  the  homogeneous  crust,  and 
hence  the  transmission  of  Lg  is  still  fairly  efficient. 
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Figure  50.  Seismic  sections  of  model  3.  which  has  a  moderate  free-surface  topography  and  a  flat  Moho.  The  scattering 
effects  of  topography  on  Lg  are  not  as  strong  as  those  due  to  large-scale  structural  discontinuities. 


62 


FDstep  0 

VEFTT  .OOOsec,  1  800  1  500 
XTIcMark=  10  km;  Z  Tic  Mark  =  10  km 
-O^E-fOO,  0^E4^ 


FDstep  400 

VERTIOXXXIaac,  1  800  1  500 
XTicMai1c=  10  km;  Z  Tic  Marks  10km 
-Oil5E+00,  O^E+00 


FDstep  800 

VERTZO^sec,  1  800  1  500 
XTIcMarks  10  km;  Z  Tic  Marks  10km 
-035E-i^,  O^E+00 


FDstep  1200 

VERT30.000sec,  1  800  1  500 
XTIcMarks  10  km;  Z  Tic  Mark  s  10km 
-031E+00,  030E^ 


FDstep  1600 

VERT40i)00sec,  1  800  1  500 
XTIcMailcs  10  km;  Z  Tic  Marks  10km 
^J29E4Q0,  0.27E400 


FDstep  2000 

VERT50.000SOC,  1  800  1  500 
XTIcMarks  10  km;  Z  Tic  Marks  10km 
-0^4E^,  O^E-^00 


FDstep  2400 

VERT60.000sec,  1  800  1  500 
XTIcMarks  10  km;  Z  Tic  Mark  s  10km 
•O^E-01,  O^E-01 


LFD  Simulation  of  Lg  Propagation:  Model  3a 

Figure  51.  The  vertical-component  snapshots  of  Lg  wave  propagation  of  model  3a,  in  which  the  Moho  is 
fully  compensated.  The  elevation  change  of  this  topographic  profile  is  relatively  small  compared  to  the 
thickness  of  the  homogeneous  crust,  and  hence  the  transmission  of  Lg  is  still  fairly  efficient. 
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Rgure  52.  Seismic  sections  of  model  3a,  which  has  a  moderate  free-surface  topography  and  a  fully  compensated 
Moho.  The  scattering  effects  of  topography  on  Lg  are  not  as  strong  as  those  due  to  large-scale  structural  discontinuities. 
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LFD  Simulation  of  Lg  Propagation:  Model  4 

Figure  53.  The  vertical-component  snapshots  of  Lg  wave  propagation  of  model  4. 
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Figure  54.  Seismic  sections  of  model  4,  which  has  a  trapezoidal  topography  and  a  fully  compensated  Moho. 
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are  not  the  same  overtone  of  the  Rayleigh  mode. 


5  801  10  801  15  801 

Model  OC,  0km  0.135E+00  0-P  Model  OC,  10km  0.159E+00  0-P  Model  OC,  40km  (-Sn)  0.123E-01  O-P 
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Figure  56.  Same  as  Figure  55  except  for  models  OC,  OCa,  OD,  ODa,  and  ODb.  All  these  models  have  embedded  sedimentary  layer(s)  and  hence  the  Lg  energy  is 
delayed.  There  is  relatively  less  energy  that  would  penetrate  the  Moho,  as  shown  on  the  waveforms  at  40  km.  The  model  ODb  generates  enhanced  Rg  waves 
(bottom  left). 
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Figure  57.  Same  as  Figure  55  except  for  models  OE,  OEa,  OEb,  OEc,  and  OEd.  Models  of  the  “E”  series  excite  very  strong  Sm  phase. 


Model  OF,  40km  (~Sn)  0.569E- 
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Figure  58.  Same  as  Figure  55  except  for  models  OF ,  OFa,  OFb,  2a,  and  4.  The  moderate  surface  and  Moho  topography  in  model  2a  affect  Lg  waveform  charac¬ 
teristics  only  mildly.  Model  OFb  excites  very  strong  Rg  phase,  very  similar  to  model  ODb  (see  bottom  left  of  Figure  56). 


4.  L  PROPAGATION  IN  ANELASTIC  MODELS 

g 

The  incorporation  of  intrinsic  attenuation  due  to  an  arbitrary  absorption  law  is  believed  to  be 
straightforward  in  frequency-domain  methods.  However,  so  far  many  codes  used  in  waveform  synthesis 
(such  as  the  reflectivity  method  and  the  wavenumber-integration  method)  can  only  handle  a  frequency- 
independent  Q.^  Bache  et  a!.  (1981)  found  that  Lg  synthetics  made  for  models  with  frequency- 
independent  Q  do  not  attenuate  with  the  proper  frequency  dependence.  If  models  are  constructed  that 
match  Lg  amplitude-distance  relationship  at  1  Hz,  then  Lg  synthetics  would  attenuate  too  rapidly  at 
higher  frequencies. 

For  the  time-domain-based  numerical  methods,  it  has  been  very  difficult  to  add  the  intrinsic 
attenuation  because  the  anelastic  stress-strain  relation  has  the  form  of  a  convolution  integral,  which  is 
intractable  in  a  numerical  computation.  Vidale  and  Helmberger  (1988)  convolved  finite-difference  syn¬ 
thetics  with  a  time-varying  operator  to  model  the  effect  of  anelastic  attenuation.  This  method  is  not  suit¬ 
able  for  media  with  spatially-varying  Q.  This  method  is  also  not  appropriate  when  there  are  significant 
wave  conversions  (P  to  S  etc. )  in  media  where  Qp  is  not  equal  to  Qs-  The  first  successful  attempt  to 
incorporate  realistic  attenuation  laws  into  time-domain  methods  was  made  by  Day  and  Minster  (1984) 
based  on  the  method  of  Fade  approximation,  which  yields  an  expansion  of  the  frequency-dependent 
viscoelastic  modulus  into  a  rational  function.  Emmerich  and  Kom  (1987)  propose  a  slightly  different 
method  based  on  the  rheological  model  of  the  generalized  Maxwell  body,  which  has  a  modulus  of  the 
desired  rational  form.  The  m^or  inconvenience  of  these  approaches  is  that  they  demand  a  continual 
storage  of  five  or  more  time  steps  of  the  wave  field,  depending  on  the  accuracy  of  the  approximation. 
The  non-causal  nature  of  some  of  these  techniques  is  another  fundamental  drawback,  at  least  conceptu¬ 
ally,  since  in  reality  the  anelastic  attenuation  of  the  Earth  should  act  in  a  causal  manner.  That  is,  the  dis¬ 
sipation  of  energy  should  occur  as  soon  as  the  seismic  wave  arrives,  and  that  the  resulting  displacement 
at  the  current  time  step  should  not  be  dependent  on  future  displacement. 

With  all  these  considerations  in  mind,  a  different  procedure  is  developed  in  this  study  to  incor¬ 
porate  the  anelastic  attenuation.  It  turns  out  that  if  we  drop  the  ambitious  attempt  of  imposing  an  arbi¬ 
trary  (that  is,  user-defined)  frequency-dependence  on  the  Q  operator,  then  it  becomes  very  easy  to 
implement  a  causal,  phase-independent  damping  operator  which  is  quite  suitable  for  the  LFD  calcula¬ 
tion.  Several  researchers  have  readily  demonstrated  that  the  performance  of  commonly  used  absorbing 
boundary  conditions  can  be  greatly  improved  if  a  viscous  damping  zone  is  added  to  the  grid  boundary 
(Ceqan  et  a!.,  1985;  Levander,  1985a).  The  damping  zone  simply  reduces  the  amplitude  in  a  pointwise 
manner.  There  is  no  reason  why  this  technique  can  not  be  exploited  to  model  the  anelastic  attenuation. 

’  A  wavenumber-integration  code  with  this  limitation  lifted  is  being  developed  at  S-Cubed  by  K.  McLaughlin. 


71 


Different  seismic  waves  are  used  as  the  initial  pulse  in  these  examples  to  demonstrate  the 
effectiveness  of  this  algorithm. 

Consider  the  simplest  isotropic  homogeneous  medium  in  which  the  amplitude  of  seismic  waves 
decays  exponentially  with  traveling  distance; 

A(f,A)  =  Ao(f)  G  =  Ao(f)  e''' ^  .  [1] 

7C  ■  f 

where  y  =  q  ,  U  is  the  group  velocity,  and  A  is  the  distance  traveled.  In  LFD  calculations,  A  is  taken 
to  be  the  distance  that  the  seismic  wave  would  travel  within  one  temporal  step  of  the  LFD  iteration,  that 

_  n  -t  •  dl 

is,  A  =  U  •  dt.  Thus  G  can  also  be  written  as  e  °  .  There  exists  a  constant  ti  such  that 

G  =  e‘’  "^  =  1 -ti  y-A.  [2] 

It  implies  that  G  =  1  -  ^  ^  ^itcI  hence  Q  =  Qq  -  f  if  we  define  Oq  to  be  If  the  damp¬ 

ing  term  G  is  a  function  of  the  grid  coordinate  only  and  invariant  with  frequency,  then  we  would  have  a 
Q  inaeasing  linearly  with  frequency.  If,  however,  a  frequency-dependent  G  is  used  at  separate  LFD 
simulations,  then  combining  the  band-limited  LFD  results  would  produce  the  solution  for  that  particular 
frequency-dependent  Q  model.  Here  the  parameter  q  is  a  function  of  y  or  G,  and  the  means  of  deter¬ 
mining  q  will  be  discussed  later.  In  practice,  however,  the  users  only  need  to  specify  a  multiplicative 
constant  G  slightly  less  than  1  for  each  grid  point.  These  localized  damping  factors  are  used  to  modify 
the  displacement  field  pointwise  at  each  iteration  step.  The  decay  rate  (y)  and  the  quality  factor  (Q)  can 
be  determined  later,  after  the  finite-difference  calculation  is  done.  A  possible  drawback  of  this  approach 
is  that  if  a  specific  frequency-dependent  Q  model  is  desired,  then  several  separate  LFD  simulations 
need  to  be  carried  out  for  each  frequency-0  pair  of  interest,  as  discussed  above.  Nevertheless,  this 
possible  shortcoming  of  this  approach  is  outweighed  by  its  simplicity.  More  importantly,  this  procedure 
presences  the  causality.  Another  characteristic  of  this  approach  is  that,  given  a  damping  factor  G  in  a  P- 

_ja 

SV  LFD  calculation,  Qq  =  ^  would  be  applicable  to  both  P  and  S  phases.  Thus,  Op  and  Os 

should  be  about  the  same. 


So  far  we  have  derived  several  necessary  conditions  for  an  anelastic  attenuation  model,  based  on 
the  desired  exponential  decay  of  seismic  amplitude.  In  the  following,  we  shall  take  a  schematic  view  of 
this  proposed  algorithm.  Consider  the  heterogeneous  acoustic  wave  equation  in  the  nondissipative 
medium; 


a^p 

df 


^^(x,  2)[ 


a^p 

ax2 


2p 

az2 


[3] 


where  P  is  the  acoustic  displacement  potential,  and  c(x,  z)  is  the  acoustic  velocity  at  the  node  (x,  z). 
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LFD  iteratively  solves  for  the  current  pressure  at  (x,  z)  using  that  of  the  past  within  a  spatially-limited 
region  surrounding  the  grid  (x.  z).  The  actual  size  of  the  temporal  and  spatial  buffers  required  depends 
on  the  order  of  the  LFD  scheme.  At  each  time  step,  the  resulting  unattenuated  pressure  P  is  then  multi¬ 
plied  by  the  damping  G  to  get  the  dissipated  P.  Once  the  whole  pressure  field  is  updated  and  dissi¬ 
pated  (grid  by  grid)  with  the  corresponding  spatial  damping  factors,  the  standard  LFD  iteration  restarts 
for  the  next  pressure  field  without  damping,  and  then  damping  factors  are  applied  again.  The  procedure 
for  the  P-SV  LFD  calculation  is  exactly  the  same.  It  should  be  obvious  that  the  extra  calculation  and 
memory  required  in  this  procedure  are  minimum.  We  can  envision  the  pointwise  damping  factor  as  a 
degenerate  digital  filter  which  has  only  one  point  of  temporal  span  (that  is,  memory).  It  can  be  regarded 
as  the  limiting  case  of  the  convolutional  integral  of  many  attenuation  operators  typically  used  in  the 
frequency-domain  approach. 

The  attenuated  pressure,  P  =  P  -  G,  can  also  be  written  as 

p  ^  p  _  n  •  .  p  [4] 

That  is,  P  can  be  obtained  by  adjusting  the  undissipated  P  a  little.  The  coefficient  of  pressure  loss, 

is  very  similar  to  the  coefficient  of  friction  term  that  Levander  (1985)  (and  Frankel  and  Wenner- 
Q 

berg,  1987)  used  in  the  telegraphy  equation.  The  difference  is  that  we  have  added  the  t]  term  here  to 
account  for  the  correct  relationship  between  Q  and  the  damping  effect.  Jih  (1996)  gives  a  more  detailed 
discussion  as  well  as  several  examples  using  different  seismic  waves  as  input  to  demonstrate  the 
effectiveness  of  this  algorithm. 

Seven  models  with  Qq  values  varying  from  60  to  500  are  tested  in  this  study.  The  attenuative  por¬ 
tion  is  100  km  long  and  30  km  thick.  In  each  case,  a  significant  fraction  of  the  energy  is  lost  due  to 
absorption.  Nevertheless,  the  chessboard  pattern  of  the  initial  Lg  wave  packet  is  retained  for  all  models. 
In  fact,  none  of  these  seven  models  exhibits  any  notable  scattering  (see  Figures  59  through  62).  Also, 
the  linear  frequency-dependence  of  the  resulting  Q  is  quite  clear,  so  is  the  frequency-independent  y 
(Table  3).  These  observations  strongly  indicate  that  the  simple  attenuation  operator  proposed  for  the 
LFD  method  is  adequate  for  operational  purposes. 

Following  the  same  quantification  procedure  as  in  the  previous  section,  amplitude  ratios  represent¬ 
ing  the  Lg  transmission  and  coupling  to  other  phases  are  computed.  Instead  of  using  the  surface  syn¬ 
thetic  to  measure  the  quality  factor  and  y,  the  amplitude  ratio  of  transmitted  Lg  waves  is  based  on  that 
at  a  depth  of  10  km,  which  we  find  to  be  a  more  stable  measure  of  the  transmitted  energy.  The  results 
are  summarized  in  Table  3  below. 
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Table  3.  Lg  Propagation  Statistics  of  7  Anelastic  Models 


Model,  e"'^^ 

T 

Qo 

Y{1  Hz) 

Q(2  Hz) 

y(2Hz) 

OQa,  0.99747 

0.175 

0.200 

0.053 

0.026 

0.014 

59 

0.0162 

139 

0.0138 

OQb,  0.99873 

0.418 

0.445 

0.049 

0.026 

0.013 

117 

0.0082 

249 

0.0077 

OQc,  0.99916 

0.590 

0.614 

0.056 

0.027 

0.014 

196 

0.0049 

398 

0.0048 

OQd,  0.99937 

0.672 

0.695 

0.059 

0.028 

0.015 

263 

0.0036 

531 

0.0036 

OQe,  0.99949 

0.718 

0.739 

0.060 

0.028 

0.017 

316 

0.0030 

635 

0.0030 

OQf,  0.99958 

0.770 

0.060 

0.028 

0.018 

396 

0.0024 

796 

OQh,  0.99968 

0.826 

0.835 

0.060 

0.029 

0.020 

530 

0.0018 

1054 

0.0018 

BK:  spectral  amplitude  ratio  at  1  Hz  of  “on-time  Lg  ”  (3.33  km/sec)  to  that  of  the  reference  model. 

T:  spectral  amplitude  ratio  at  1  Hz  of  transmitted  Lg  (recorded  at  a  depth  of  10  km)  to  that  of  the  reference  model. 

Sn/Lg;  spectral  amplitude  ratio  at  1  Hz  of  converted  to  the  reference  Lg  . 
far-P/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  teleseism'c  P  to  the  reference  Lg  . 
far-S/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  teleseismic  S  to  the  reference  Lg  . 

As  Q  increases,  the  transmission  effidency  (T)  also  increases,  as  expected.  However,  the  coupling 
of  Lg  into  S„  or  teleseismic  phases  does  not  seem  to  be  affected.  Redudng  the  thickness  of  the 
attenuating  body  does  not  seem  to  change  this  observation  (Figure  62).  This  is  different  from  the  case 
of  Rg  in  which  Jih  (1995,  1996)  reported  that  the  undissipated  Rg  energy  could  actually  be  converted 
into  S  or  Lg  waves.  The  reason  for  this  difference  is  not  dear  as  yet.  A  possible  explanation  is  that  it 
has  to  do  with  the  fundamental  difference  between  Rg  and  the  higher  modes.  The  partide  motion  of  Rg 
is  very  localized  in  that  the  energy  causing  retrograde  (or  prograde)  rolling  at  greater  depths  does  not 
interact  with  the  shallow  portion  of  the  crust.  When  the  rolling  at  shallow  depths  is  dissipated,  the 
“waveguide”  is  actually  changed  and  the  undissipated  energy  at  greater  depth  can  no  longer  propagate 
in  the  same  mode  as  .  As  a  result,  H  is  converted  into  a  pure  shear  or  Lg  -like  wave  (see  Jih,  1995). 
On  the  other  hand,  Lg  waves  are  the  interference  of  planar  S  waves  multiply  reflected  at  the  free  sur¬ 
face  and  the  Moho.  This  mode  of  propagation  is  more  robust  in  that  even  if  part  of  the  energy  on  the 
wavefront  is  dissipated  in  the  upper  crust,  the  remaining  part  of  the  wavefront  might  still  be  able  to  pro¬ 
pagate  more  or  less  along  the  original  direction.  Lg  energy  in  the  deeper  crust  will  also  be  dissipated 
when  it  eventually  reaches  the  upper  crust.  Likewise,  the  Lg  energy  in  the  shallow  portion  of  the  crust 
will  be  dissipated  first  and  then  eventually  enter  the  deeper  crust  and  so  forth.  Thus  the  shallow 
attenuating  layer  will  affect  the  whole  Lg  wavefield  in  a  relatively  more  uniform  manner  without  breaking 
up  the  mode  (eigenfunction).  As  a  result,  there  will  be  no  notable  enhancement  of  S„  or  teleseismic 
coupling. 
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Figure  59.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  an  anelastic  mode,  OQd,  with 
Oo  {Lg )  =  263  for  Lg  waves.  The  anelastic  portion  of  the  model  is  100  km  long  and  30  km  thick.  In  this 
case,  the  effect  of  the  intrinsic  attenuation  is  to  reduce  the  amplitude  through  dissipation.  No  significant 
scattering  is  observed. 


LFD  Simulation  of  Lg  Propagation:  Model  OQd, 122 
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Synthetics  at  0km  Depth.  Model  0C3d,l22 

VERT,  same  scale,  decimated  by  5 


0.117E+00 


Synthetics  at  0km  Depth,  Model  0Qd,122 

HORI,  same  scale,  decimated  by  5 


Figure  60.  Seismic  sections  of  model  0Qd,l22  which  has  an  attenuating  body  100  km  long,  30  km  thick,  with  Qo(f^ )  = 
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amplitude  but  does  not  cause  any  mode  conversion. 
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LFD  Simulation  of  Lg  Propagation:  Model  0Qd,62 

Figure  62.  Similar  to  Rgure  59  except  a  15km-thick  attenuative  layer.  For  the  Lg  phase,  the  effect  of 
intrinsic  attenuation  is  to  reduce  the  amplitude  through  dissipation.  Changing  the  volume  size  (that  is,  the 
thickness)  of  the  attenuating  body  does  not  seem  to  cause  more  mode  conversion.  This  is  very  different 
from  the  case  for  Rg . 
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5.  L  PROPAGATION  IN  MODELS  WITH 
RANDOM  VELOCITY  VARIATION 

In  an  attempt  to  model  the  Lg  blockage  observed  in  the  west  Pyrenees  Range,  Chazalon  et  al. 
(1993)  and  Gibson  and  Campillo  (1994)  suggest  that  scattering  by  small-scale  features  within  the  lower 
crust  Is  the  primary  reason  for  the  observed  blockage.  In  this  section,  several  heterogeneous  models 
with  various  RMS  random  velocity  fluctuations  are  tested  to  quantify  the  scattering  effects  of  small-scale 
heterogeneity  on  Lg .  Figure  63  shows  the  vertical-component  snapshots  of  Lg  wave  propagation  in  a 
single-layer  model  with  a  heterogeneous  portion  100  km  long  and  30  km  thick  embedded  in  it.  The  self¬ 
similar  random  heterogeneity  in  this  sandwiched  portion  has  an  RMS  velocity  fluctuation  of  10  percent 
and  a  scale  length  of  1  km.  The  random  medium  is  generated  with  the  procedure  described  In  Frankel 
and  Clayton  (1986). 

As  the  RMS  velocity  variation  increases  to,  say  20  percent,  the  scattered  field  becomes  very  com¬ 
plicated  {cf.  Rgures  63  and  65).  Contrary  to  the  case  with  intrinsic  attenuation  alone,  the  small-scale 
random  heterogeneity  causes  many  kinds  of  coupling  (scattering)  (Figure  65).  Table  4  below  summar¬ 
izes  all  the  results.  The  “blocking”  parameter  (BK)  was  computed  at  two  different  depths,  0  and  10  km, 
and  the  latter  is  used  in  determining  Qq  and  y.  Both  the  transmission  of  Lg  and  the  coupling  of  other 
phases  are  directly  affected  by  the  change  in  the  RMS  velodty  variation.  As  the  RMS  velocity  variation 
increases,  the  amplitude  ratio  of  transmitted  Lg  to  the  reference  Lg  wave  drops  significantly.  The  scat¬ 
tered  energy  enhances  the  coupling  of  Sp  as  well  as  those  of  teleseismic  phases.  It  is  interesting  to 
note  that  Lg  -to-S„  coupling  is  not  a  simple  linear  function  of  RMS  velocity  variation,  unlike  the  coupling 
of  teleseismic  phases  or  Lg  transmission.  Figure  68  and  Table  4  indicate  that  below  10%,  the  Lg  -to-Sn 
coupling  does  exhibit  a  highly  linear  trend  with  the  increasing  complexity  of  the  medium.  The  coupling 
drops  as  the  RMS  variation  goes  beyond  10%.  Intuitively  this  could  be  because  the  scattered  wavefield 
would  become  dominant  and  the  isotropic  point  scatterers  would  radiate  the  energy  in  all  directions. 

Another  observation  is  that  as  the  RMS  variation  increases,  the  delay  of  the  centroid  of  the  Lg 
wave  packet  on  the  surface  synthetics  (see  Rgure  67)  becomes  very  obvious.  This  is  similar  to  the  sto¬ 
chastic  dispersion”  phenomenon  McLaughlin  and  Anderson  (1987)  reported  for  P  waves.  In  theory,  the 
apparent  attenuation  and  the  dispersion  should  form  a  Hilbert  transform  pair,  according  to  the  Kramers- 
Kronig  relation  (Jacobson,  1987).  This  can  be  tested  if  the  frequency  band  Is  broad  enough. 
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Table  4.  Lg  Propagation  Statistics  of  Random  Models 


Model 

BK(0) 

BK(10) 

SnfLg 

far-P/Lg 

far-S/Lg 

Q(1  Hz) 

y(1  Hz) 

OY,  2% 

1.020 

0.978 

0.065 

0.040 

0.042 

3934 

0.0002 

OY,  3% 

1.032 

0.952 

0.075 

0.048 

0.055 

1878 

0.0005 

OY,  4% 

1.044 

0.917 

0.086 

0.057 

0.069 

1088 

0.0009 

OY,  5% 

1.059 

0.873 

0.097 

0.065 

0.083 

707 

0.0014 

OY,  6% 

1.077 

0.822 

0.108 

0.073 

0.096 

494 

0.0019 

OY.  7% 

1.095 

0.766 

0.118 

0.081 

0.108 

365 

0.0026 

OY,  8% 

1.114 

0.706 

0.126 

0.089 

0.118 

280 

0.0034 

OY.  9% 

1.132 

0.644 

0.133 

0.097 

0.127 

222 

0.0043 

OY.  10% 

1.146 

0.583 

0.138 

0.106 

0.134 

181 

0.0053 

OY,  11% 

1.154 

0.524 

0.141 

0.115 

0.140 

150 

0.0064 

OY,  12% 

1.154 

0.471 

0.143 

0.126 

0.144 

128 

0.0075 

OY.  13% 

1.141 

0.423 

0.143 

0.136 

0.147 

111 

0.0086 

OY,  14% 

1.115 

0.381 

0.142 

0.145 

0.150 

98 

0.0098 

OY.  15% 

1.075 

0.344 

0.138 

0.154 

0.155 

88 

0.0108 

OY,  16% 

1.020 

0.313 

0.133 

0.161 

0.161 

81 

0.0118 

OY,  17% 

0.949 

0.285 

0.126 

0.167 

0.167 

75 

0.0127 

OY,  18% 

0.859 

0.260 

0.118 

0.171 

0.173 

70 

0.0136 

OY,  19% 

0.749 

0.235 

0.110 

0.174 

0.177 

66 

0.0146 

OY.  20% 

0.625 

0.210 

0.103 

0.174 

0.181 

61 

0.0158 

BK(0):  spectral  amplitude  ratio  at  1  Hz  of  “on-time  Lg  ”  (3.33  km/sec)  to  that  of  the  reference  model. 

BK(1 0):  spectral  amplitude  ratio  at  1  Hz  of  transmitted  Lg  (recorded  at  a  depth  of  1 0  km)  to  that  of  the  reference  model. 

Sn/Lg:  spectral  amplitude  ratio  at  1  Hz  of  converted  Sp  to  the  reference  Lg  . 

far-P/Lg:  spectral  ampRtude  ratio  at  1  Hz  of  converted  teleseismic  P  to  the  reference  Lg  . 

far-S/Lg;  spectral  amplitude  ratio  at  1  Hz  of  converted  teleseismic  S  to  the  reference  Lg  . 


Both  anelasticity  and  small-scale  random  heterogeneity  can  contribute  to  the  so-called  “Lg  block¬ 
age”.  Note  that  in  practice  the  Lg  blockage  occurs  whenever  the  peak  Lg  amplitude  is  below  the  signal 
level  of  other  phases  surrounding  Lg .  It  certainly  does  not  require  the  signal  level  in  the  predicted  Lg 
window  to  be  far  below  the  noise  level.  In  our  calculations,  a  velocity  variation  of  8  percent  or  higher 
leads  to  an  equivalent  spatial  Q  of  270  or  less,  which  would  be  sufficient  to  reduce  the  Lg  amplitude  by 
30  percent  or  more  for  every  100  km  it  traverses  (see  the  parameter  T  in  Table  3).  The  classical  exam¬ 
ple  of  “Lg  blockage”  along  the  path  from  the  Novaya  Zemlya  test  site  to  northern  Scandinavia  has  a 
Qo{Lg)  about  252,  as  discussed  before.  Based  on  Table  4,  this  level  of  attenuation  is  approximately 
equivalent  to  a  model  with  a  velocity  variation  of  8  percent  alone.  Stable  shield  regions  are  reported  to 
have  an  RMS  velocity  variation  less  than  5  percent.  For  such  paths,  the  corresponding  spatial  Qo  value 
would  be  700  or  larger. 
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Some  comparisons  of  various  conversion  mechanisms  may  be  interesting.  For  a  mountainous 
model.  Lg  -to-S  is  twice  as  strong  as  Lg  -to-P  (Table  2).  (Here  P  and  S  refer  to  the  teleseismic  P  and  S, 
respectively.)  For  the  “A"  and  “B”  series  of  models,  as  well  as  the  random  media,  these  two  mechan¬ 
isms  are  about  equally  strong  (see  Tables  2  and  4).  Anelastic  models  and  “C”-type  models  have  ^  Lg- 
to-S  weaker  than  the  Lg  -to-P  (Table  3),  whereas  models  of  the  “D”,  “E”,  and  “F”  series  show  the  oppo¬ 
site  trend  (Table  2). 
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LFD  Simulation  of  Lg  Propagation:  Model  0Y,10% 

Figure  63.  The  vertical-component  snapshots  of  Lg  wave  propagation  in  a  single-layer  model  with  a 
heterogeneous  portion  100  km  long  and  30  km  thick  embedded  in  it.  The  self-similar  random  heterogeneity 
in  this  sandwiched  portion  has  an  RMS  velocity  fluctuation  of  10  percent  and  a  scale  length  of  1  km. 
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Rgure  64.  Seismic  sections  of  the  model  OY.10%.  This  model  has  10  percent  velocity  variation  in  an  embedded 
heterogeneous  volume  which  is  1 00  km  long  and  30  km  thick. 
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LFD  Simulation  of  Lg  Propagation:  Modei  0Y,20% 

Figure  65.  The  verticai-component  snapshots  of  Lg  wave  propagation  in  a  single-layer  model  embedded 
with  a  heterogeneous  portion  100  km  long  and  30  km  thick  in  it.  The  self-similar  random  heterogeneity  in 
this  sandwiched  portion  has  an  RMS  velocity  fluctuation  of  20%  and  a  scale  length  of  1  km.  The  strong 
scattering  generates  a  very  complicated  wavefield  (as  compared  to  Figure  63).  Much  of  the  scattered 
energy  would  contribute  to  coda  waves. 
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Rgure  66.  Seismic  sections  of  the  model  0Y,20%.  This  modei  has  20%  velocity  variation  in  an  embedded  heterogene¬ 
ous  volume  which  is  100  km  long  and  30  km  thick.  This  model  generates  very  long  coda,  as  expected.  Also  note  that 
the  centroid  of  the  Lg  wave  packet  is  delayed. 
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Figure  67.  Synthetics  recorded  at  3  depths  (0km.  10km,  and  40km)  after  crossing  a  heterogeneous  zone  100  km  wide  and  30  km  thick.  As  the  velocity  fluctuation 
increases  from  2  percent  (top)  to  20  percent  (bottom),  the  Lg  amplitude  decreases  (as  expected).  The  centroid  of  the  Lg  wavetrain  also  gets  delayed  somewhat. 
The  Lg  -to-Sn  coupling  decays  once  the  RMS  variation  is  greater  than  13-15  percent,  presumably  because  most  of  the  energy  has  already  been  lost  to  scattering  in 
all  directions.  A  variation  of  2  percent  leads  to  very  little  scattering  and  hence  the  synthetics  are  very  similar  to  those  of  the  reference  model  (see  Figure  61). 
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pling  to  teleseismic  phases  become  more  prominent. 


6.  REMARKS  ON  OBSERVATIONS  OF  Lg  BLOCKAGE 
AND  A  MODEL-BASED  CORRECTION  PROCEDURE 

Baumgardt  (1991)  conducted  a  thorough  observational  study  of  Lg  blockage  for  Soviet  explosions 
recorded  at  regional  or  far-regional  distances.  He  documented  several  examples  of  partial  and  complete 
blockage  and  carefully  related  the  possible  cause  of  each  case  to  geological  and  geophysical  charac¬ 
teristics  along  the  path.  Baumgardt  found  that  in  every  instance  of  partial  or  complete  blockage,  the 
great-circle  path  of  Lg  between  the  source  and  the  receiver  aosses  contained  sedimentary  basins  and 
adjacent  topographic  discontinuities.  Lg  blockage  does  not  seem  to  relate  to  variations  in  total  crustal 
thickness.  Blockages  result  from  basins  where  there  are  unusually  thick  sediments,  usually  in  excess  of 
10  km,  where  the  basin  is  fully  contained  and  the  sediments  pinch  out  at  the  edges  of  the  basin,  and/or 
where  there  are  topographic  inhomogeneities  at  the  edges  of  the  basin.  He  observed  less  blockage  for 
events  occurring  inside  sedimentary  basins  themselves,  where  the  basins  are  not  particularly  thicK  or 
where  the  sediments  thin  very  gradually  with  distance.  Table  4  below  summarizes  observations  of  Lg 
blockage  and  weakening,  along  with  the  possible  mechanism  that  has  been  proposed. 


Table  5  Observed  Lg  Blockage/Weakening  and  Proposed  Mechanism(s) 


Mechanism  Proposed 

Authorfs) 

Region  Studied 

Dramatic  thinning  of  waveguide 

Ewing  et  al.  (1957) 

continent-ocean  margin 

Pure  oceanic  paths 

Knopoff  et  al.  (1979) 

Atlantic  &  Pacific  Oceans 

Missing  granitic  layer 

Piwinskii  (1981) 

Caspian  Sea 

Crustal  thickening 

Variations  in  crustal  thickness  (+Q?) 
Variations  in  crustal  thickness 

Ruzaikin  et  al.  (1 977) 
Gregersen  (1984) 

Kennett  et  al.  (1985) 

Himalayan  Belt 

North  Sea  Graben 

North  Sea  Graben 

Anelastic  attenuation  in  sediments 

Mitchell  &  Hwang  (1987) 

U.S. 

Scattering  by  small-scale  heterogeneity 
Scattering  by  small-scale  heterogeneity 

Campillo  et  al.  (1 993) 
Chazalon  et  al.  (1993) 

SW  Alpine  Range 

Western  Pyrenean  Range 

Scattering  from  tectonic  boundaries 
Scattering  from  tectonic  boundaries 
Scattering  from  tectonic  boundaries 

Kadinsky-Cade  et  al.  (1981) 
Ni  &  Barazangi  (1 983) 
Baumgardt  (1985, 1990) 

Turkish  &  Iranian  Plateaus 

Tibetan  Plateau 

Ural  Mountains 

Contained  thick  sedimenteuy  basin 
Contained  thick  sediments  (+low  Q?) 

Baumgardt  (1991) 

Baumgardt  (1991) 

Barents  Shelf 

North  Caspian  Depression 
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The  importance  of  understanding  crustal  waveguide  effects  on  Lg  waves  has  motivated  many 
modeling  studies  of  blockage  in  recent  years  (Kennett,  1986a:  Maupin,  1989;  Regan  and  Harkrider, 
1989;  Chazalon  et  al.,  1993;  Gibson  and  Campillo,  1994;  and  Cao  and  Muirhead,  1993).  Some  of 
these  authors  suggest  that  the  observed  blockage  is  not  straightforward  to  simulate,  and  that  the  unmo¬ 
deled  small-scale  random  heterogeneity  is  required  to  account  for  the  blockage,  on  top  of  the  large-scale 
structural  (that  is,  geometrical)  variation  of  the  waveguide.  Although  either  low  Q  or  large  velocity  varia¬ 
tion  can  be  invoked  to  explain  the  Lg  blockage  of  many  paths,  we  might  be  able  to  narrow  down  the 
possible  mechanisms  with  some  additional  diagnostic  information.  For  instance,  if  a  pronounced  Lg  -to- 
Sn  conversion  is  observed  along  with  the  Lg  blockage,  then  it  is  likely  that  the  blocking  mechanism  is 
not  due  to  low  Q. 

One  the  m^or  challenges  in  monitoring  a  Comprehensive  Test  Ban  Treaty  [CTBT]  is  discriminating 
small  events  at  regional  distances.  Various  forms  of  compressional/shear  (P/S)  ratios  have  been  pro¬ 
posed  and  extensively  tested  as  candidate  discriminants.  In  tackling  such  problems,  a  better  amplitude 
(or,  equivalently,  magnitude)  scale  of  the  regional  phase  of  interest,  with  respective  path  effects  carefully 
accounted  for,  is  definitely  useful  or  sometimes  even  necessary.  Along  this  line  of  thought,  Jih  et  al. 
(1995)  (and  also  Jih  and  Lynnes.  1993)  suggest  a  simple,  convenient  magnitude  scale  for  Lg : 

=  4.0272  -  Bias  +  logA(A)  +  ^log(A(km))  +  ^log[sin(^^^  .  [51 


where  Y=^7g.  Q(0  =  Qo  ' 

A  is  the  epicerrtral  distance  in  km  and  A(A)  is  the  observed  raw  Lg  amplitude  measured  in  the  time 
domain  in  pm  [microns]  at  the  epicerrtral  distance  of  A  km.  Or,  equivalently. 


A(km) 


K  =  m  ■  ■  <A(kn.))=  ■  t)1 
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^-KA-IOkm) 


[Sal 


The  “Bias”  term  in  Eqs.  (5)  and  (5a)  is  meant  to  account  for  the  different  Lg  excitation  (relative  to  mb).  It 
is  set  to  zero  for  the  Eastern  U.S.  Thus  a  seismic  source  with  1-sec  Lg  amplitude  ofllOpmatlOkm 
(extrapolated)  epicerrtral  distance  would  correspond  to  a  nttg  of  4.0272  +  2.0414  +  0.3333  -  1.4019  + 
0.0000  =  5.000,  which  has  been  suggested  to  be  appropriate  for  both  Eastern  North  America  and  Semi- 
palatinsk.  That  is  to  say,  a  seismic  source  in  these  two  regions  with  mb  5.0  would  have  a  approxi¬ 
mately  the  same.  Jih  et  al.  (1995)  suggest  use  of  0.34  and  0.26  m.u.  for  the  “Bias”,  respectively,  for 
Pahute  Mesa  and  Novaya  Zemlya  explosions.  For  earthquakes  in  the  Iranian  Plateau,  Nuttli  (1980)  sug¬ 
gested  a  bias  value  of  0.39  m.u.,  which  includes  both  the  bias  due  to  regional  tectonics  as  well  the 
Lg  -mb  bias  due  to  source  type. 
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If  the  regionalized  y  map  is  available  (such  as  those  in  Singh  and  Hemmann,  1983;  Jih  and 
Lynnes,  1993;  and  many  others)  a  path-specific  y  for  an  arbitrary  source-station  pair  can  simply  be  com¬ 
puted  as  the  weighted  sum  of  the  yk’s  of  the  subregions  that  the  ray  path  traverses: 

K  Ak 

y=2:yk-x-  16] 

However,  it  is  dear  from  Eq.  (2)  that  an  erroneous  path  y  would  yield  a  bias  that  increases  with  the 
distance.  Furthermore,  this  error  is  independent  of  the  actual  source  size  or  the  quality  of  the  raw  ampli¬ 
tude  measured  at  the  recording  station.  It  is  the  bias  solely  due  to  inaccurate  calibration  of  the  propaga¬ 
tion  effect.  Jih  and  Lynnes  (1993)  gave  several  examples  showing  how  Lg  amplitude  measurements 
could  be  severely  biased  if  a  poor  estimate  of  the  path  y  is  used  or  if  the  path  y  is  not  used  at  all. 
Baumgardt  (1995)  recently  demonstrated  that  by  incorporating  a  distance-correction  term,  which  is 
essentially  the  ^  in  Eq.  (5a),  Into  his  Pg  iLg  ratio  computation,  a  much  better  discrimination  result 
was  obtained.  In  fact,  this  is  exactly  the  same  reason  why  Nuttli’s  (1986, 1988)  absolute  yield  estimates 
of  underground  nudear  explosions  using  patfi-corrected  are  good. 

The  seismic  discrimination  problem  is  actually  intertwined  with  the  source-size  determination  prob¬ 
lem,  despite  the  misperception  that  the  estimation  of  source  size  in  general  (and  the  estimation  of  explo¬ 
sive  yield  in  particular)  may  no  longer  be  as  Important  in  the  CTBT  context.  P/S-ratlo-based  discrimina¬ 
tion  is  a  procedure  to  discern  differences  in  energy  partitioning  among  phases  and  among  frequency 
bands.  Miscellaneous  amplitude  measurements  needed  to  identify  the  event  are  actually  reflecting 
different  parts  of  the  focal  sphere  and/or  frequency  contents  of  the  seismic  source  itself.  These  parame¬ 
ters  are  often  obscured  by  the  path  effects.  Thus,  a  better  path-calibration  procedure  would  always 
benefit  the  discrimination,  for  just  the  same  reason  it  would  benefit  the  source-size  determination.^ 

In  a  situation  where  some  geometrical  blockage  occurs  along  the  path,  that  is,  there  is  a  major 
structural  variation  along  the  path,  it  may  be  possible  to  further  refine  Eq.  (5).  As  an  example,  the  r(A- 
10km)/ln(10)  term  for  the  path  from  Novaya  Zemlya  to  KEV  can  be  broken  down  into  yi  (Ai-10km)/ln(10) 
-I-  y2-A2/ln(10).  where  y2  =  0.0052  (that  is,  Qo=185;  see  the  model  OD  in  Table  2),  A2  is  the  lateral  span 
of  the  Terrigenous  Sediments  (Baumgardt,  1991),  and  y,  is  the  average  Lg  attenuation  along  the  great- 
circle  path  excluding  the  sedimentary  segment.  Obviously  it  is  necessary  to  acquire  relevant  geological 
and  geophysical  information  along  the  path  before  applying  this  model-based  correction. 


^  See  discussion  in  Dainty  (1995)  and  Blandiord  et  al.  (1992). 
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7.  DISCUSSION  AND  CONCLUSIONS 


An  accurate  prediction  of  the  regional  phases  in  areas  of  high  proliferation  concern  requires  a 
decent  understanding  of  the  attenuation  and  scattering  mechanisms  along  the  propagation  paths.  Syn¬ 
thetic  data  have  been  of  great  benefit  in  gaining  physical  insight  into  discriminants,  and  may  be  the  only 
means  of  evaluating  evasion  scenarios.  Synthetic  seismograms  are  particularly  useful  for  regions  where 
earthquake  or  explosion  data  are  not  available.  There  is  no  doubt  that  both  geometrical  blockage  (due  to 
large-scale  lateral  structural  variation),  Intrinsic  attenuation,  and  small-scale  random  heterogeneity  in  the 
crust  affect  Lg  propagation.  But  it  is  often  very  difficult  to  separate  these  effects  observalionally.  Again, 
numerical  modeling  could  be  the  only  means  to  quantitatively  identify  the  respective  roles  played  by 
these  factors. 

In  this  study,  the  two-dimensional  LFD  method  has  been  used  to  model  the  propagation  and 
scattering  of  Lg  waves  in  a  suite  of  crustal  models,  including  rough  interfaces,  rugged  free-surface 
topography,  and  embedded  thick  sedimentary  layers.  Two  non-standard  techniques  have  been 
specifically  developed  for  this  study:  [1]  a  pure  Lg  wave  packet  (for  triggering  LFD  calculations),  and  [2] 
a  crude,  causal  attenuation  operator.  The  pure  Lg  wave  packet  is  injected  into  a  stratified  portion  of  the 
grid  as  the  reference  initial  condition  to  trigger  all  LFD  calculations.  This  is  a  natural  extension  of  Jih’s 
(1993b)  work  on  Rg  to  Lg  problems.  In  addition  to  those  commonly  recognized  advantages  of  LFD 
method,  another  obvious  advantage  of  this  approach  is  that  the  effect  due  to  different  types  of  hetero¬ 
geneity  on  the  seismic  phase  of  interest  can  be  isolated  and  evaluated  easily.  The  other  tool  tested  in 
this  study  is  a  new  attenuation  operator  suitable  for  time-domain  numerical  calculations.  As  a  crude 
approximation,  it  applies  to  only  a  narrow  band  within  which  it  gives  a  nearly  frequency-independent  y. 
as  expected.  The  most  notable  advantage  Is  its  simplicity  to  implement.  Furthermore,  it  is  a  causal 
operator  which,  at  least  conceptually.  Is  more  realistic.  Overall  it  appears  to  perform  reasonably  well. 
M^or  observations  made  with  LFD  calculations  are  summarized  in  the  following: 

Early  L  :  The  Moho  uplift,  which  is  common  in  models  of  the  “A”  and  “F’  series,  breaks  down  the 
Lg  waveguide  and  causes  Lg  -to-S„  coupling  at  the  ascending  interlace.  When  the  converted  S„  waves 
encounter  the  descending  Moho  of  these  models,  some  energy  penetrates  the  interface  again,  either  as 
a  leg  of  Sn  or  as  converted  Lg .  The  Moho  uplift  in  “F”-series  models  is  not  as  dramatic  as  In  A  - 
series  models.  Nevertheless,  both  the  Lg  -to-Sn  and  possibly  S„  -to-L^  conversions  are  still  apparent 
The  early  Lg  observed  for  paths  from  Novaya  Zemlya  to  Scandinavia  can  be  explained  with  this 
mechanism.  The  fact  that  Lg  can  be  excited  at  continental  margins  from  oceanic  S„  conversion  has 
long  been  observed  and  reported  (for  example,  Isacks  and  Stephens,  1975). 
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Late  Lg  and  longer  coda  waves:  The  sedimentary  layer,  which  is  present  in  all  “C”-series  models 
and  some  of  “D”  and  “F”  models,  is  the  apparent  cause  of  late  Lg  amVals.  For  models  ODb,  OFb  and 
OFc,  the  uppermost  soft  layer  causes  strong  reverberations  as  well  as  a  very  strong  Lg  -to-Rg  conver¬ 
sion.  By  seismic  reciprocity,  there  should  be  some  structures  that  would  cause  strong  Rg  -to-Lg  conver¬ 
sion.  Part  of  this  issue  has  been  addressed  in  Jih  (1995).  Regarding  the  coda  decay  rate,  intrinsic 
attenuation  and  scattering  (particularly,  that  by  random  heterogeneity)  appear  to  have  different  effects. 
Increased  scattering  attenuation  causes  more  energy  to  be  distributed  into  the  coda  from  the  direct 
phase  with  increasing  time.  On  the  other  hand,  inaeased  intrinsic  attenuation  entails  a  greater  loss  of 
elastic  energy  with  time  from  both  the  coda  and  the  direct  wave.  This  observation  can  be  verified  by 
comparing  synthetic  seismograms  in  Figures  61  and  67,  and  is  in  agreement  with  Frankel  and  Wenner- 
berg  (1987). 

excitation  and  S^/L^  ratios:  Except  for  models  00,  2,  2a,  3,  3a,  OC,  OCa,  and  OCb,  all  other 
models  excite  significant  S„  waves  through  the  Lg  -to-S„  conversion.  A  thick  sedimentary  layer  by  rtsetf 
does  not  necessarily  cause  strong  Lg  -io-Sn  conversion.  Models  of  the  “C”-series  are  such  examples. 
However,  if  the  thick  sediment  layer  is  contained,  such  as  in  models  of  “D”  series,  then  some  Sn  waves 
can  be  generated  at  the  pinched  interface.  Models  with  irregular  Moho  topography  generate  Sn  with  a 
different  mechanism,  as  discussed  above.  Most  models  excite  converted  S„  waves  in  the  forward 
direction.  Hoever,  models  OA  and  OB  also  excite  some  Sn  in  the  backward  direction.  The  extremely 
large  Sn  /Lg  ratio  associated  with  models  OE,  OEa,  and  OEc  (all  with  a  thickened  cnist)  offers  an  excel¬ 
lent  demonstration  for  the  observation  Ruzaikin  et  al.  (1977)  made. 

Topographic  relief:  Models  with  a  mild  free-surface  topography  alone  (for  example,  models  2  and 
3)  generate  a  relatively  weak  Lg  -to-Sn  conversion  comparable  to  that  of  “C”  models.  Comparing  the 
synthetics  in  Figures  55  and  58,  it  appears  that  the  transmitted  Lg  wave  trains  of  models  with  mild 
mountainous  topography  are  very  similar  to  those  of  the  reference  model  00.  In  fact,  these  mountainous 
models  and  models  with  weak  random  heterogeneity  as  well  as  anelastic  models  would  stand  out 
among  all  models  in  that  the  checkerboard-like  interference  pattern  of  the  original  Lg  wave  train  is 
(somewhat)  retained  (see  Rgures  45,  49,  59,  and  62).  This  is  very  different  from  the  case  of  Rg  for 
which  the  rugged  free-surface  topography  is  shown  to  be  a  strong  scatterer  (Jih,  1993b,  1995). 

Lg  to  Pg  conversion:  All  models  of  the  “F’  series  exhibit  significant  Lg  -io-Pg  conversion.  The 
peak  horizontal  amplitude  on  trace  No.  25  ranges  from  10  to  15  percent  of  that  of  the  original  Lg  wave 
train.  To  a  lesser  extent,  models  of  “C”  and  “D”  series  as  well  as  models  OA  and  OBa  also  generate  the 
same  conversion. 

Coupling  to  teleseismic  phases:  For  mountainous  models,  Lg  -to-S  is  twice  as  strong  as  Lg  -to- 
P.  For  the  “A”  and  “B”  series  of  models,  as  well  as  the  random  media,  these  two  mechanisms  are 
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about  equally  strong.  For  anelastic  models  and  “C”-type  models,  Lg  -to-S  is  weaker  than  the  Lg  -to-P, 
whereas  models  of  the  “D”,  “E”,  and  “P’  series  show  the  opposite  trend  (Table  2). 

Lg  blockage:  Both  abrupt  changes  in  the  Moho  topography  and  a  thick  contained  sedimentary 
layer  can  cause  strong  Lg  -\o-S„  conversions.  The  physical  mechanisms  underlying  the  conversion  are 
different  for  these  two  structures,  however.  The  early  Lg  observed  for  paths  from  Novaya  Zemlya  island 
to  Scandinavia  is  illustrated  with  LFD  calculations.  Alluvial  basin  structures  cause  strong  reverberations 
as  well  as  very  strong  Lg  -to-Rg  conversions.  Valley  structures  with  large  depth  to  width  ratios  and 
large  velocity  contrast  with  the  substrate  support  very  large  amplitude  lateral  resonance  modes.  If  we 
invoke  the  principle  of  seismic  reciprocity,  Rg  -to-Lg  conversion  could  be  an  important  mechanism  of  Lg 
excitation  for  certain  structures.  Recent  observational  study  by  Patton  and  Taylor  (1995)  suggests  that 
this  is  the  case,  at  least  for  Yucca  Rat  explosions. 

Other  than  major  structural  discontinuities  in  the  crustal  waveguide,  both  anelastidty  and  small- 
scale  random  heterogeneity  can  also  contribute  to  the  so-called  “Lg  blockage”.  For  instance,  an  RMS 
velocity  variation  of  8  percent  in  the  whole  crust  is  equivalent  to  a  Qo{Lg )  of  about  270,  which  would  be 
sufficient  to  reduce  the  peak  amplitude  of  IHz  Lg  waves  by  30  percent  for  every  100  km  they  traverse. 
The  classical  example  of  an  Lg  -blocking  path  between  Novaya  Zemlya  and  northern  Scandinavia  has  a 
Qo  {Lg )  about  252  and  could  be  modeled  with  an  RMS  velocity  variation  of  8  percent  alone. 

Unlike  Rg  waves,  Lg  prop^ation  is  not  affected  by  a  moderate  free-surface  topography  alone  as 
much  as  other  types  of  heterogeneity.  Topographic  relief  has  been  shown  to  be  a  very  strong  scatterer 
for  the  Rg  phase  (Jih,  1993b,  1995).  This  is  because  Rg  energy  is  confined  in  the  uppermost  layer. 
The  energy  of  Lg  is  evenly  distributed  in  the  crust,  and  the  surface  topography  variations  are  generally 
small  relative  to  the  crustal  thickness.  Thus  it  may  be  expected  that  the  topographic  variations  alone 
should  have  relatively  weaker  influence  on  Lg  propagation,  as  compared  to  Rg .  This  is  illustrated  with 
models  2  and  3  of  this  study.  However,  the  topography  may  reflect  the  crustal  thickness  through  isos- 
tasy.  Zhang  and  Lay  (1994)  suggest  that  S„  iLg  ratios  may  be  correlated  to  some  roughness  measures 
of  the  paths,  based  on  a  limited  data  set.  LFD  appears  to  be  the  ideal  tool  to  test  the  existence  of  any 
correlation  between  some  roughness  measure  and  the  propagation  characteristics  of  regional  phases, 
as  demonstrated  with  the  very  few  mountainous  models  tested  in  this  study.  Similar  statistical  analysis 
should  also  be  conducted  to  examine  the  effects  of  water  column,  random  Moho  topography,  and  irregu¬ 
lar  intra-crustal  interfaces,  as  well  as  the  combination  of  all  these  factors. 
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